NPS  ARCHIVE 

1969 

CYR,  J. 


INVESTIGATION  OF  THE  SENSITIVITY 

OF  A  PATROLLING  SUBMARINE'S  CAPABILITY 

TO  GAIN  A  DETECTION  AS  A  FUNCTION  OF 

INCREASING  SURFACE  SHIP  SPEEDS 


by 


Joseph  Henry  Cyr 


TZ- 


j- 


United   States 
Naval  Postgraduate  School 


THESIS 


INVESTIGATION  OF  THE  SENSITIVITY  OF  A 

PATROLLING  SUBMARINE • S 

CAPABILITY  TO 

GAIN  A  DETECTION  AS  P 

i  FUNCTION  OF 

INCREASING  SURFACE 

SHIP  SPEEDS 

by 

Joseph  Henry 

Cyr 

and 

Leonard  Bento 

Santos 

October  1969 


Tku>  document  ka&  bzen  appmovtd  ^on.  pub  tic  kz- 
Iza&z  and  haJLz;  itb  dUtxibution  ti>  unUmitzd. 


\\  \  ' }  7  7 


Investigation  of  the  Sensitivity  of  a 
Patrolling  Submarine's  Capability  to 
Gain  a  Detection  as  a  Function  of 
Increasing  Surface  Ship  Speeds 


by 


Joseph  Henry  Cyr 
Lieutenant  Commander,  United  States  Navy 
B.S.,  Purdue  University,  1960 

and 

Leonard  Bento  Santos 

Lieutenant  Commander,  United  States  Navy 

B.A.,  Brown  University,  1959 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  OPERATIONS  RESEARCH 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
October  1969 


ABSTRACT 


Numerical  integration  and  Monte  Carlo  techniques  are 
used  in  the  development  of  several  models  in  order  to 
determine  the  effect  on  probability  of  random  detection 
of  a  merchant  ship  using  speeds  up  to  90  knots  by  a  10 
knot  submarine  patrolling  a  back-and-forth  barrier.   A 
definite  range  law  for  detection  is  assumed.   Individual 
encounter  models  are  developed  for  ship  tracks  which  cross 
the  midpoint  of  the  submarine  patrol  line  at  various  angles. 
The  models  are  extended  to  include  the  assumption  of  a 
normal-distribution  of  crossing  points.   Computer  programs 
of  the  models,  written  in  the  FORTRAN  IV  language,  are  in- 
cluded.  The  results  are  applied  in  a  numerical  example. 

It  is  concluded  that  while  increases  in  ship  speeds  do 
result  in  a  substantial  decrease  in  probability  of  detection 
by  a  submarine  in  the  case  of  a  single  barrier  transit,  a 
speed  advantage  alone  when  applied  to  a  typical  transit  of 
the  North  Atlantic  will  not  appreciably  decrease  the  overall 
detection  probability. 
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I.   INTRODUCTION 

A.   BACKGROUND 

1.   High-Speed  Ship  Development 

New  principles  and  concepts  in  ship  hull  design 
have  made  possible  the  development  of  ships  capable  of 
speeds  heretofore  considered  impossible.   In  particular, 
the  advent  of  surface-effect  ship  (air  cushion  and  captured 
bubble),  hydrofoil,  and  catamaran  hull  designs,  coupled 
with  advances  in  nuclear  power,  gas  turbine,  and  water- jet 
propulsion  systems,  will  have  revolutionary  effects  on 
naval  strategy  and  tactics,  and  world  shipping  patterns. 
That  serious  consideration  to  high-speed  merchantmen  is 
being  made  is  evidenced  by  the  fact  that  the  U.S.  Navy  and 
the  Maritime  Administration  share  funding  in  the  Joint 
Surface-Effect  Ship  Program  which  was  established  in  1966 
to  determine  the  feasibility  of  4,000  to  5,000-ton  surface- 
effect  ships  capable  of  speeds  in  excess  of  80  knots.   Bell 
Aerosystems,  Inc.,  and  Aerojet-General  Corp.  have  each 
received  $1.5  million  as  a  start  on  a  3-year  program  to 
build  an  experimental  100-ton  surface-effect  ship  [1] .   The 
fact  that  high-speed  merchant  ships  which  can  enjoy  large 
speed  advantages  over  submarines  are  indeed  in  the  offing 
indicates  a  need  for  reexamination  of  the  abilities  of 
future  submarines  to  detect  and  kill  these  ships.   This 
thesis  makes  a  start  by  examining  the  effect  of  ship- 
submarine  speed  ratios  on  detection  probability 
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resulting  from  random  encounters  by  submarines  on  fixed 
patrols . 

a.   Literature  Search 

As  might  be  expected,  the  literature  dealing 
with  high-speed  ships  is  relatively  sparse;  those  papers 
that  do  make  reference  to  ships  of  the  future  concern 
themselves  with  the  hydrofoil  small-ship  types  [2],  or  with 
the  so-called  super  merchant  ships  which  have  such  features 
as  large  (500 , 000-ton)  displacements,  and  underwater 
tankers  of  speeds  up  to  40  knots  [3].   Various  Navy  studies 
in  the  past  which  have  attended  to  the  question  of  ship 
speeds  versus  vulnerability  have  been  limited  in  scope  to 
surface  ship  speed  ranges  up  to  only  35  knots  [4,  5].   One 
study  in  particular  concluded  that  "...When  ships  are  sailed 
independently ...  (and)  when  cost  is  considered  there  is  only 
a  marginal  gain  in  going  from  20  to  24  knots...."  /I/.   The 
cost  consideration  at  that  time  (1955)  involved  studies  of 
conventional  hulls.   An  earlier  study  [6]  showed  a  decrease 
in  percentage  of  independent  merchant  ships  sunk  from  10-12 
per  cent  when  ship  speeds  were  approximately  equal  to 
submerged  submarine  speeds  (4-8  knots)  to  1-2  per  cent  when 
ship  speeds  were  double  that  of  the  U-boat  submerged  speeds. 
When  hull-mounted  sonars  on  merchant  ships  were  considered 
[7]  ,  maximum  optimal  sustained  speeds  for  independent 
vessels  were  in  the  range  of  20  to  24  knots. 

In  summary,  neither  the  search  of  available 
studies  [8],  operations  research  papers,  nor  computer 
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searches  for  documentation  [9]  at  the  Defense  Documentation 
Center,  Washington,  D.C.,  and  at  the  U.S.  Naval  Postgraduate 
School,  Monterey,  California  have  revealed  any  literature 
regarding  high-speed  ships,  and  in  particular  the  effect  of 
merchant  ship  speed  increases  beyond  35  knots  on  vulner- 
ability to  submarines.   Indeed,  all  facets  associated  with 
these  high-speed  ships  should  be  examined,  e.g.,  noise 
radiation  in  water,  the  submarine  fire  control  problem, 
effect  of  torpedoes,  etcetera. 

2 .   Vulnerability  of  Convoys  and  Independent  Ships 

As  was  clearly  evidenced  during  World  Wars  I  and  II, 
convoys  and  independent  merchants  were  quite  susceptible  to 
attack  by  submarines.   Allied  shipping  losses  from  the 
German  U-boat  campaign  during  the  first  five  years  of  World 
War  II  provide  sufficient  evidence  of  the  extreme  vulner- 
ability of  merchant  shipping  [10,  11].   Considering  the 
improvements  in  both  submarine  and  submarine  weapon  system 
capabilities  since  this  era,  the  shipping  losses  previously 
experienced  can  be  regarded  as  the  minimum  losses  that  may 
be  expected  in  a  future  war  [10].   Furthermore,  it  is  a 
foregone  conclusion  that  in  the  opening  stages  of  a  future 
general  war  escort  numbers  will  be  insufficient  to  provide 
adequate  convoy  protection  for  all  merchant  shipping  [12, 
13],   It  can  be  assumed,  then,  that  the  initial  phases  of 
a  future  war  will  see  an  aggressive  enemy  submarine  force 
presented  with  many  unescorted  targets  which  can  be  attacked 
with  impunity  [10]. 
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To  counter  the  possibility  of  sustaining  heavy 
merchant  losses  at  the  outset  of  a  future  war,  passive 
defense  measures  such  as  installation  of  sonar  on  merchant 
ships,  routing  to  minimize  losses,  convoying,  and  increas- 
ing merchant  ship  speeds  have  been  considered  [12].   This 
thesis  considers  the  passive  defense  measure  of  increased 
merchant  ship  speed. 

3 .   Future  Capabilities  of  Submarines 

The  Type  XXI  submarine  produced  by  the  Germans  in 
1944  influenced  basic  submarine  design  until  the  appearance 
of  the  nuclear-powered  submarine  in  1954  [14].   Since  the 
introduction  of  the  nuclear-powered  submarine,  great  tech- 
nological strides  have  been  realized  in  the  areas  of  sub- 
marine hull  construction,  propulsion,  sonar,  torpedo 
weapons  systems,  and  fire  control  systems.   These  advances 
have  contributed  significantly  to  the  threat  to  merchant 
shipping  posed  by  submarine  forces  of  potential  enemies. 
The  modern  nuclear-powered  submarine  can  equal  or  exceed 
the  speed  of  its  surface-ship  target  or  pursuer,  and  can 
remain  submerged  for  days.   Together  with  this  endurance 
capability,  the  nuclear-powered  submarine  is  able  to  main- 
tain a  high  sustained  submerged  speed  in  a  relatively 
quiet  environment.   "...The  limiting  lines  of  approach  to 
convoys  and  task  forces  steadily  widen  with  increasing 
submarine  speed  until  there  is  no  sector  safe  from  attack 
by  a  30-knot  submarine...."  /2/    This  quote  was  taken  from 
a  publication  dated  April  30,  1962;  since  this  time  a  great 
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deal  of  research  and  development  in  the  area  of  surface- 
effect  ships  has  taken  place.   The  Joint  Surface-Effect 
Ship  Program  was  initiated  to  consider  large  surface- 
effect  ships  capable  of  high  speeds.   It  is  conceivable 
that  in  the  forseeable  future  surface  ships  with  speeds 
far  exceeding  the  capabilities  of  the  nuclear-powered 
submarine  will  take  their  place  in  the  surface  ship  com- 
munity.  Submarine  speeds  are  not  expected  to  exceed  their 
present  maximum  for  some  time  to  come,  thus  permitting 
the  high-speed  merchant  ship  to  enjoy  a  comfortable  ad- 
vantage in  speed. 

4 .   Safety  of  Independents 

Heretofore  convoys  have  been  recommended  as  the 
best  alternative  for  moving  merchant  vessels  across  the 
oceans  when  a  submarine  threat  exists  [15],   This  is  due 
primarily  to  the  fact  that  for  slow  moving  (less  than  15 
knots)  merchants,  safety  lies  in  numbers.   The  measure  of 
effectiveness  for  analysis  of  merchant  shipping  has  been 
percentage  lost  of  the  total  number  of  ships  in  the 
convoy  [16].   Studies  completed  during  World  War  II  indi- 
cated that  approximately  the  same  number  of  ships  were  lost, 
per  attack,  from  large  convoys  as  from  small  convoys.   Hence 
the  conclusion  was  to  not  only  continue  convoying,  but  to 
increase  the  size  of  the  convoys  as  well. 

With  the  advent  of  high-speed  surface-effect  ships, 
the  problem  of  independent  sailing  versus  convoying  should 
again  be  considered.   In  a  study  conducted  by  the 
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Antisubmarine  Operations  Research  Group  Tenth  Fleet  in 
December  194  3,  the  advantages  and  disadvantages  that  might 
result  from  independent  sailings  were  discussed:  "....Fast 
merchant  vessels,  with  the  exception  of  those  carrying  high 
priority  cargo,  have  been  able  to  continue  sailing  inde- 
pendently along  the  coast  (Atlantic)  much  of  the  time."  /3/ 
Here  is  a  direct  reference  to  a  merchant  ship  speed  advan- 
tage over  the  submarine. 

Other  advantages  of  independent  sailings  are  of 
two  kinds:   first,  the  independent  merchant  can  operate  more 
efficiently  and  hence  deliver  more  cargo  in  a  given  time; 
second,  the  reduced  requirements  for  escorts  releases  these 
surface,  air,  and  submerged  craft  for  offensive  roles  [11]. 
5 .   Probabilities  of  Encounter 

The  problem  of  random  encounter  between  two  units, 
sometimes  referred  to  as  the  searcher  and  the  target  [17, 
18,  19,  20],  is  examined  in  most  of  the  classical  operations 
research  literature,  Operations  Evaluation  Group  reports 
and  studies  (which  constitute  a  large  part  of  the  so-called 
"classical  literature"),  and  more  recently  by  search  theor- 
ists and  game  theorists. 

Probabilities  of  encounters,  or  detections  resulting 
from  applications  of  search  plans,  are  discussed  in  much  of 
the  literature,  ranging  from  Kimball  and  Morse  [16]  and 
Koopman  [18],  through  the  summary  reports  of  the  Operations 
Evaluation  Group  [14]  and  applications  studies  [2]  to  the 
modern  search  theorists,  such  as  Pollock  [19,  20]. 
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Data  concerning  encounters  versus  ship  position 
displacement  from  reported  submarine  contacts  [21]  and 
frequency  of  encounters  versus  submarine  position  predic- 
tions [22]  have  been  gathered  and  analyzed  during  the 
Second  World  War.   Applications  of  game  theory  for  the 
selection  of  convoy  routes  which  would  minimize  ship- 
submarine  encounters  were  proposed  in  a  study  by  Danskin 
of  the  Operations  Evaluation  Group  in  1954  [15]  ,  which 
later  appeared  as  a  paper  in  Operations  Research  [23] .   The 
Random  Encounter  with  no  involved  search  effort  has  been 
studied  by  several  authors;  Koopman  [17,  18]  considered  the 
random  encounter  of  an  observer  progressing  at  a  constant 
velocity  within  a  uniform  distribution  of  targets  of  con- 
stant speed,  as  well  as  within  a  circular-normal  distribu- 
tion of  targets.   The  problem  of  random  encounter  between 
a  surface  force  progressing  at  constant  velocity  with 
another  unit  travelling  at  a  fixed  speed,  the  position  and 
heading  of  which  are  independently  and  normally  distributed, 
has  been  examined  by  Dobbie  in  1945  [24]. 

Literature  concerning  the  effect  of  speed  on  the 
random  encounter  is  limited  almost  exclusively  to  government 
studies  and  reports.   Kittel,  in  1944,  analyzed  data  con- 
cerning sinkings  and  sightings  by  submarines  to  give  an 
early  indication  of  the  effect  of  speed  on  the  safety  of 
ships  [6].   The  results  were  discussed  by  Sternhell  and 
Thorndike  in  1946  [25],  and  by  Winston  in  1955  [4].   These 
results,  in  turn,  were  summarized  by  the  Operations 
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Evaluation  Group  in  1957  [12].   More  recently,  the  general 
problem  of  the  effect  of  speed  on  the  vulnerability  of 
independent  merchant  vessels  was  examined  by  Neufer  in 
1961  [4].   Indeed,  Neufer  appears  to  have  been  the  first 
to  attempt  to  determine  the  explicit  survival  probability 
of  a  merchant  vessel  in  transit  from  one  port  to  another  as 
a  function  of  its  speed. 

Finally,  the  probability  of  encounter  by  a  unit  on 
a  barrier  patrol  has  been  examined  by  at  least  two  authors. 
Koopman,  in  his  Search  and  Screening,  examined  the  problem 
using  the  underlying  assumption  of  a  uniform  distribution 
of  points  of  intersection  of  the  target's  track  with  the 
back-and- forth  barrier  patrol  line.  /4/   In  particular, 
Koopman  showed  that,  with  a  definite  range  law,  the 
probability  of  detection  by  the  back-and- forth  barrier 
patrol  is  given  by 


1  -  [A  -  (/l/62  +  1  -  l)/2]  /A  (A  +  1), 


for  6>  l/(2/A  (A  +  1)'  ) 


1,     for  6  <  l/(2/A  (A  +  1)'  ) 

where 

A  =  ratio  of  the  barrier  patrol  line  length  to 

sweep  width 
6  =  ratio  of  target  speed  to  barrier  unit  patrol 

speed 
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Further,  Koopman  showed  that,  for  all  values  of  A,  the 
back-and- forth  barrier  patrol  is  preferred  to  the  crossover 
("figure  eight")  type  patrol  for  all  values  of  6  £.  1.  /5/ 
Van  Train  [26]  used  a  war  gaming  method  in  a  limited  way  to 
determine  the  probability  of  detection  of  a  transitor 
intercepting  a  multi-unit,  single  or  multi-line  fixed 
barrier,  in  which  each  submarine  was  constrained 

1)  to  remain  in  the  center  of  an  assigned  zone,  and 

2)  to  conduct  intercepts  using  headings  parallel  to  the 
barrier  line.   The  submarine's  course,  and  position  on  a 
"starting  line"  were  assumed  to  have  been  normally- 
distributed. 

B.   THE  PROBLEM 

The  basic  encounter  problem  is  that  of  determining  the 
probability  of  detection  of  a  surface  ship  by  a  submarine 
on  an  assigned  back-and-forth  barrier  patrol,  such  as  that 
discussed  by  Koopman  [18] .   The  following  assumptions  were 
made  for  the  models: 

1.   Assumptions  for  the  Basic  Random  Encounter  Model 

a.  The  submarine  patrols  a  geographically-fixed 
back-and-forth  barrier  centered  on  a  known  surface  ship 
track  or  convoy  route . 

b.  The  submarine  patrols  at  a  constant  speed,  v. 

c.  A  definite  range  law  for  detection  is  used. 

d.  The  submarine  patrol  length,  L,  is  four  times 
the  detection  radius. 
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e.  The  submarine  position  at  any  time,  x,  is  from 
a  uniform  distribution. 

f.  The  point  of  intersection  of  the  ship's  track 
and  the  patrol  line  (or  extension  of  the  patrol  line)  is 
from  a  normal  distribution,  with  standard  deviation  of  l/6th 
of  the  length  of  the  submarine  patrol  line. 

C.   METHODS  OF  SOLUTION 

Two  methods  of  solution  to  the  models  were  devised;  the 
first  was  a  numerical  integration  technique,  i.e.,  simula- 
tion by  systematic  sampling.   The  submarine's  initial 
position  was  successively  incremented  by  500  yards  (0.25 
nautical  mile)  and  the  minimum  lateral  range  (CPA)  between 
the  submarine  and  the  surface  vessel  was  computed  by 
numerical  vector  analysis.   The  basic  model  considered  960 
initial  submarine  positions  for  each  speed  ratio  considered; 
twelve  ship-submarine  speed  ratios ,  ranging  from  1  to  9  were 
examined.   The  basic  models  also  examined  nine  angles  of 
intersection  of  tracks,  ranging  from  90  down  to  30  degrees. 
The  extended  model  considered  ten  specific  points  of  inter- 
section and  eight  speed  ratios. 

The  second  method  may  be  classed  as  a  Monte  Carlo  method, 
at  least  within  Brown's  definition:   "...any  procedure  which 
involves  the  use  of  statistical  sampling  techniques  to 
approximate  the  solution  of  a  ...  physical  problem. ..."  /6/ 
The  Monte  Carlo  method  utilized  virtually  the  same  computer 
program  subroutines  as  the  numerical  integration  approach, 
but  selected  the  submarine's  initial  positions  from  a 
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uniform  pseudorandom  number  (URN)  generator  [27]  ,  and 
selected  track  intersection  points  using  the  Box-Mueller 
method  of  normal  random  number  generation  [28].   For  this 
technique,  the  basic  model  took  1,000  samples  for  each  of 
nine  speed  ratios  per  crossing  angle,  and  for  17  crossing 
angles.   The  extended  model  took  10,000  samples  for  each 
of  nine  speed  ratios.   The  total  number  of  vector 
( "maneuvering-board" )  solutions,  then,  exceeded  500,000; 
the  usefulness  of  the  computer  in  this  analysis  is  obvious. 
1 .   The  Random  Encounter  Model 

Consider  a  submarine  patrolling  a  barrier  of  length 
L  miles  at  a  constant  speed,  v.   Assume  that  the  position 
of  the  submarine  on  its  patrol  line  at  any  time,  x ,  is  a 
uniformly  distributed  random  variable,  Y,  .   Without  loss  of 
generality,  then,  the  problem  can  be  considered  to  begin  at 
time  t  when  a  surface  vessel  on  a  course  a,  projected  to 
intersect  the  submarine  patrol  line  at  its  center,  is  a 
distance  x  /cos  a  from  the  point  of  intersection.   Define 
the  distance  x   to  be  equal  to  the  range  of  the  submarine 
detection  (assuming  a  definite  range  law) ,  R,  which  in  the 
initial  case  was  considered  to  have  been  30  miles.   Refer 
to  Figure  1  and  the  definitions  summarized  below: 

->  ? 

u  =  u  cos  a  i  +  u  sin  a  j    =  ship  velocity  ; 

-*■       1 

v  =  ±  v  J  =  submarine  velocities; 

-»■->■-> 

w,  =  v,  -  u  =  relative  velocity  of  the 

submarine  when  it  is  on  the 
up leg; 
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—  (x   ,y  +p  ) 


(xp'yb} 


Submarine 

Patrol 

Line 

■  VV 


Figure  1.   Vector  Diagram  of  the 
Encounter  (Relative  Space) 
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w„  =  v„-  u  =  relative  velocity  of  the 

submarine  when  it  is  on  the 
down leg; 

t,  =  (y   +  L  -  Y,)/v      =  time  remaining  for  the 
"  submarine  to  be  on  its 

initial  leg  before  course 
reversal; 


t,  =  (Y,  -  y  ) /v  =  time  remaining  for  the  sub- 

p  marine  to  be  on  its  initial 

downleg  before  course 
reversal; 


T  =  L/v  =  length  of  time  required  for 

the  submarine  to  travel 
entire  length  of  its  patrol 
line; 

±      ?      *  ... 

R,  =  x  1  +  Y  j  =  initial  position  vector  of 


1     p      b- 


the  submarine  from  the  ship; 


R  =  R-.    +   wn  t  =  position  of  the  submarine 

mil  i.    •  a.   i    •  j. 

when  it  makes  its  course 

reversal. 
2 .   Algorithm  to  Determine  Probability  of  Encounter 

(1)  Select  the  initial  value  of  Y,  such  that  Y,  =  y  . 

b  b    J  p 

(2)  Specify  the  initial  direction  of  submarine  velocity; 
call  it  v. 

(3)  Compute  the  values  of  t,  t,  ,  R,  ,  and  R;  set  w  =  w,  . 

(4)  Compute  the  dot  product  R-, -w  =  z: 

(a)  If  z  >  0,  CPA  =  R, ;  go  to  step  (5); 

(b)  If  z  <  0,  compute  the  dot  product  R  *w  =  z1; 

1)  If  z1  <  0,  CPA  =  R  ; 

—  '  n 

2)  If  z'  >  0,  the  CPA  lies  between  R,  and  R  ,    and 
must  be  determined; 

a)  set  a  =  1/2; 

b)  compute  y  =  (aw)/|w|; 

c)  compute  the  dot  product  (R,  +  y) «w  =  z"; 
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1.)  if  z"  <0,  add  1/2  to  a,  and  go  to 
step  2)b)  above; 

2.)  if  z"  >  0,  then  CPA  =  R,  +  y;  go  to 
step  (5) . 

(5)  (a)   If  | CPA |  <  R,  a  detection  occurred;  go  to  step 
(1)  if  the  number  of  samples  taken  is  greater  than  N, 

otherwise  go  to  step  (6). 

(b)   If  |CPA|<  R,  there  was  no  detection  on  this  leg: 

->■->•       -> 

1)  if  w  =  w,  ,  set  w  =  w~  ; 

-»■     -»■  -v     -> 

2)  if  w  =  w~ ,  set  w  =  w,  ; 

3)  go  to  step  (4) . 

(6)  P  (Encounter)  =  #  detections/  N. 

The  general  flow  diagrams  describing  the  logic  for  the 
models  are  contained  in  Appendix  B.   Note  that  there  are 
two  models,  namely,  the  Basic  Model,  in  which  the  ship  track 
crosses  the  midpoint  of  the  submarine  patrol  line  at  various 
angles,  and  the  Extended  Model,  in  which  the  crossing  angle 
is  fixed  at  90  degrees,  but  the  crossing  position  is  varied 
along  the  track  (and  its  extension) .   Appendix  C  contains 
the  actual  computer  instructions  used  to  generate  the  data 
analyzed  herein. 

D.   DESCRIPTION  OF  THE  COMPUTER  SYSTEM 

The  computer  system  used  in  the  development  of  this 
thesis  was  the  high-performance  IBM  System/360  Model  67, 
which  is  located  at  Ingersoll  Hall,  U.S.  Naval  Postgraduate 
School,  Monterey,  California.   This  sophisticated  data 
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processing  system  employs  a  FORTRAN  G  compiler,  along  with 
the  standard  software  and  hardware  associated  with  the  360 
system. 

Programming  was  accomplished  using  the  standard  FORTRAN 
IV  language  as  proposed  by  the  American  Standards  Associa- 
tion X3.4.3  FORTRAN  Working  Group  [29]. 

Approximate  execution  times  for  the  four  programs  in- 
cluded in  this  thesis  were  as  follows: 

Basic  Model:   45  minutes; 

Extended  Model:    30  minutes; 

Random  Sampling,  Midpoint:    110  minutes; 

Random  Sampling,  Normally-distributed  crossing: 
46  minutes. 
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II.   THE  RANDOM  ENCOUNTER  -  ANALYSIS 
BY  NUMERICAL  INTEGRATION  - 


A.   RANDOM  ENCOUNTER  MODEL  (FIXED  CROSSING  POINT) 
1 .   Description  of  the  Model 

The  algorithm  described  in  section  l.B.  was  used, 
except  that  the  values  of  Y,  ,  and  the  initial  direction  of 
the  submarine  were  determined  using  numerical  integration 
techniques  as  follows: 

The  submarine's  initial  position  for  the  first 
sample  was  established  at  the  bottom  of  the  patrol  line, 
with  the  initial  direction  "up."   The  relative  motion  prob- 
lem was  solved  numerically  to  determine  if  a  detection 
occurred  for  this  configuration.   For  each  subsequent 
sample,  the  submarine's  initial  position  was  incremented 
by  500  yards  (l/4t^1  of  a  nautical  mile)  ,  until  each  500 
yard  increment  in  the  initial  direction  had  been  sampled, 
up  to  the  initial  position  lying  on  the  patrol  line  extrem- 
ity.  The  submarine  direction  was  then  reversed  and  the 
procedure  continued  as  before  except  that  the  succeeding 
initial  positions  resulted  from  500-yard  increments  in  the 
new  direction.   The  number  of  samples,  then,  was  8  x  (length 
of  the  submarine  patrol  line  in  nautical  miles).   The  basic 
parameters  of  the  model  from  which  most  of  the  data  in 
Appendices  C  and  D  were  obtained  were 

submarine  patrol  line  length:   120  miles 

submarine  patrol  speed:   10  knots 

definite  range  law  detection  range:   30  miles 
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These  parameters  were  used  in  both  the  basic  and  extended 
models  for  the  simulation  by  numerical  integration,  and  for 
the  Monte  Carlo  simulations. 
2 .   Results 

The  data  resulting  from  the  simulation  by  numerical 
integration  are  contained  in  Appendix  A;  these  data  have 
been  plotted  in  Figure  2  and  Figure  3.   Figure  2a  shows  how 
the  detection  probability  varied  with  the  ship-submarine 
speed  ratio  for  a  ship  crossing  the  patrol  center  at  various 
angles.   The  points  at  which  detection  probability  was 
reduced  by  9  0%  and  95%  have  been  connected  by  the  two  curves 
labelled  "90%"  and  "95%,"  respectively.*   These  curves  are 
shown  in  Figure  2b.   Figure  3  is  an  expansion  in  scale  of 
Figure  2b,  showing  greater  detail.   Figure  4  displays  a  plot 
of  the  Probability  of  Detection  vs.  Track  Crossing  Angle  for 
a  midpoint  crossing.   The  behavior  of  P (D)  as  a  function  of 
a  is  of  little  practical  interest  for  the  case  where  a 
submarine  commander  has  freedom  to  select  a  patrol  line 
orientation  with  respect  to  a  convoy  route;  a  patrol 


it 

The  results  plotted  in  Figure  2  seems  to  indicate  that 

a  barrier  unit's  best  tactic  would  be  to  make  the  patrol  line 
cross  the  ship's  track  at  a  shallow  angle.   This  would 
indeed  be  the  case  if  the  ship's  track  were  known  with  cer- 
tainty, as  was  assumed  in  the  basic  model  discussed  here. 
The  more  realistic  assumption  of  uncertainty  in  the  know- 
ledge of  the  ship's  track  makes  the  choice  of  a  shallow 
angle  undesirable  from  the  viewpoint  of  the  barrier  unit. 
The  design  of  an  "optimal"  barrier  patrol  line,  the  crossing 
angle  of  which  depends  on,  say,  the  probability  distribu- 
tion of  a  convoy  route  displacement,  appears  to  be  an 
interesting  problem,  but  is  not  considered  herein. 
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Submarine  Speed:  10  kts . 
Patrol  Length:  120  miles 
a  =  Crossing  angle  in  degrees 
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Figure  3.   Probability  of  Random  Encounter  vs 
Ship-Submarine  Speed  Ratio  (Midpoint 
Crossing)  Enlarged  Scale 
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Figure  4.   Probability  of  Random  Encounter  vs 
Track  Crossing  Angle  (Midpoint  Crossing) 
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perpendicular  to  the  expected  surface  ship  track  would 
normally  be  selected.   Figure  4  could  be  useful  if  the 
submarine  were  restricted  for  some  reason  to  orient  its 
patrol  line  at  an  angle  other  than  90  degrees  to  the  convoy 
route.   The  effect  of  relaxing  the  assumption  that  the 
surface  ship  track  crosses  the  patrol  line  at  the  midpoint 
of  the  latter  is  examined  in  the  next  section.   The  assump- 
tion of  a  submarine  patrol  line  length  of  120  miles  (four 
times  the  detection  range)  agrees  with  Neufer's  "total 
attack  width"  which  is  derived  in  his  paper.  /I / 

B.   EXTENSION  OF  THE  MODEL  TO  VARIABLE  CROSSING  POINTS 
1 .   Description  of  the  Model 

The  basic  model  contains  the  expression 

y   =xtana  -L/2 

to  insure  that  the  ship's  track  will  bisect  the  submarine's 

patrol  line.   It  did  not  seem  unreasonable  to  assume  that  a 

submarine  normally  attempts  to  center  its  patrol  on  the 

expected  convoy  route;  that  is,  the  patrol  is  centered  on 

an  expected  value.   A  normal  distribution  of  tracks,  each 

perpendicular  to  the  submarine  patrol  line,  was  assumed  for 

the  extension  of  the  model,  so  that  Y   is  a  normally- 

P  * 

distributed  random  variable.   In  terms  of  the  coordinate 
system  used  herein, 


E(Y  )  =  -L/2 
P 
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The  logic  used  in  the  computer  program  is  illustrated  in 
Appendix  B.   The  results  of  the  computer  calculations  are 
plotted  in  Figure  5. 

2 .   Analysis  of  Results 

The  flattening  of  the  curves  for  the  speed  ratios 
p  =  1  and  p  =  9  suggested  the  possibility  of  two  simple 
limiting  shapes,  namely  the  shapes  resulting  when  p  =  0 
and  p  =  °°.   As  p  +  0,  i.e.,  when  the  ship's  speed  u  is  such 
that  u<<v,  it  can  be  seen  that  the  probability  of  detection 
will  approach  unity.   More  precisely,  for  a  given  detection 
radius,  submarine  speed,  and  submarine  patrol  line  length, 
there  exists  some  value  of  ship  speed  such  that  the  sub- 
marine will  always  have  an  opportunity  to  close  within 
detection  range.   Let 

L     =  length  of  the  submarine  patrol  line 

y     =  y-coordinate  of  the  bottom  limit  of  the  sub- 
marine  patrol  line 

u     =  ship  speed 

v     =  submarine  speed 

p     =  ship-submarine  speed  ratio 

r     =  detection  radius  (definite  range  law) 

P(D)   =  probability  of  detection 

Then,  for  all  u,  such  that 

2r/u  <  2L/v   , 

or,  for  all  p*,  such  that 

r/L  <  p*   , 
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Figure  5.   Probability  of  Detection  vs.  Crossing 
Position  (Perpendicular  Tracks) 
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lim   P(D)  =  1 


if  ly    <  L  +  r 

For  the  case  of  p  ■+  °° ,    i.e.,    submarine  speed  approaching 
zero,  the  probability  of  detection  is  simply  the  probability 
that  the  submarine's  position  (which  has  been  assumed  to  be 
uniformly-distributed  along  the  patrol  line)  falls  within 
distance  r  of  the  ship's  track.   Refer  to  Figure  6  below. 


^ 


P(D)  =  (r  +  Ay)/L 


P(D)  =  (r  -  Ay)/L 


P(D)  =  2r/L 


L 


Figure  6.   Probability  of  Detection  vs.  Crossing 
Position  for  Zero  Submarine  Speed 


It  can  be  seen  from  Figure  6  that  if  the  ship's  track 
crosses  the  patrol  line  at  a  distance  from  the  patrol  line 
end-points  greater  than  or  equal  to  r,  then  the  detection 
probability  is  simply 

2r/L 
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If  the  ship's  track  crosses  the  patrol  line  at  a  distance 
less  than  r,  say  Ay,  from  a  patrol  line  end-point,  then  the 
detection  probability  is  given  by 

(r  +  Ay)/L 

Finally,  if  the  ship's  track  does  not  cross  the  patrol  line, 
but  passes  within  detection  range  r,  say  Ay,  of  one  of  its 
end-points,  the  detection  probability  is  given  by 

(r  -  Ay)/L 


Using  the  convention  that  the  coordinate  system  origin  lies 
on  the  ship's  position  at  the  instant  that  the  ship  becomes 
"detectable,"  the  above  argument  may  be  formulated  as 
follows : 


P(D)  = 


2r/L, 


ye  (-L  +  r,  -r) 

ir 


r  -  min ( | y  I , | L  +  y  | ) sgn [ (y  ) (y   +  L) ]  , 


(1) 


y  e{[-r,r]U  [-L  -  r,  -L  +  r]  } 
y  e{[r,<»)U  (-00/  -L  -  r]} 


The  limiting  case  result,  as  presented  herein,  appears  at 
first  glance  to  disagree  with  the  limiting  case  result  ob- 
tained by  Koopman;  Koopman ' s  result  for  a  zero-speed 
barrier  unit  is  P (D)  =  0.333...  /8/,  and  the  result  from 
the  geometric  integration  for  p  =  °°  is  0.4  88.   This  latter 
value  is  a  consequence  of  the  assumption  of  normality;  if 
Koopman ' s  underlying  assumption,  namely  a  uniform 
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distribution  of  crossing  points,  were  to  be  applied  to  this 
model,  the  limiting  probability  of  detection  would  be 
(area  under  curve  for  p  =  °°)x  (density  for  uniform  distri- 
bution) or, 

(0.5)  (60)    (0.5)  (60)     60      --- 

180  180     "  180  '  u--3-3-3-"   ' 

which  is  identically  Koopman ' s  result.   The  model  expressed 
by  equations  (1)  is,  therefore,  consistent  with  Koopman. 
A  simple  FORTRAN  computer  program  was  written  (Appendix  F) 
to  compute  the  detection  probability  using  equations  (1) ; 
the  results  are  plotted  in  Figure  7.   Figure  7  also  contains 
a  plot  of  the  computed  results  for  comparison  with  the 
limiting  cases.   The  agreement  with  equations  (1)  is  evident. 
The  behavior  of  Probability  of  Detection  vs.  Crossing 
Position  in  the  region  of  ship-submarine  speed  ratios 
greater  than  three  is  unexpected  (one  would  expect  the 
curves  to  remain  uninodal) .   At  values  of  p  such  that  P>4, 
the  probability  of  detection  resulting  from  a  patrol,  the 
center  of  which  was  offset  from  the  expected  track  crossing 
by  about  l/6th  of  the  length  of  the  patrol  line,  resulted  in 
a  higher  detection  probability  than  when  the  patrol  was 
centered  on  the  track  crossing  point.   There  data  indicate 
that  a  submarine  may  effectively  reduce  the  speed  ratio  by 
properly  positioning  itself. 
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Figure  7.   Probability  of  Detection  vs.  Track 
Crossing  Point  and  Normal  Probability 
Density  of  Track  Crossing  Position 
(Perpendicular  Tracks) 
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TABLE  I 

COMPARISON  OF  PEAKS  ON  PROBABILITY  CURVES 
FOR  SHIP-SUBMARINE  SPEED  RATIOS  OF  5  AND  9 


SPEED 
RATIO 


P (D)  FOR  CENTERED 
PATROL  LINE 


P(D)  FOR  OFF-CENTERED 
PATROL  LINE  (20  MILES) 


.5104 
.5021 


.5604 
.5292 


Table  I,  extracted  from  the  data,  reveals  that  the  prob- 
ability of  detection  of  a  ship  with  a  nine-to-one  speed 
advantage  over  the  patrolling  submarine  by  a  submarine  with 
a  20-mile  offset  patrol  line  is  greater  than  the  prob- 
ability of  detection  of  a  ship  with  a  five-to-one  speed 
advantage  by  the  same  submarine  with  its  patrol  line 
centered  on  the  crossing  track.   The  table  shows  that  with 
a  20-mile  offset,  the  submarine  would  enjoy  a  higher 
probability  of  detection  than  it  would  for  the  case  of  a 
target  with  p  =  5  and  no  offset.   Figure  8  employs  a  con- 
tracted horizontal  scale  which  serves  to  emphasize  the  nodal 
behavior  discussed  above.   The  conditional  probability, 

P (Detection  I Y   =  y) , 
1  p   2 

was  determined  by  numerical  analysis,  and  is  plotted  in 
Figure  7  for  six  values  of  p.   The  track  crossing  point  was 
assumed  to  be  normally-distributed,  that  is,  N(0,400).   Now 
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Figure  8.   Probability  of  Detection  vs.  Track  Crossing 
Position  (Perpendicular  Tracks) 


42 


P  (Detection)  =  P (R  =  r)  , 

where  re  [0/00),  and  R  is  the  random  variable  representing 
the  distance  between  the  submarine  and  the  surface  ship.   In 
the  problem  at  hand,  r  =  30  miles,  and  since  a  definite 
range  law  was  assumed,  a  detection  occurs  when  R  =  30,  hence 

P (Detection)  =  P (D) 

is  a  mass  function.   From  Baye ' s  formula, 


f   y  (r,y)  =  f  (r|Y   =  y)f    (y)  .  (2) 

'  p  P       P 


Integrating  over  the  range  of  y, 


/   fpv  (r,y)dy  =  /   f(r|Y„  =  y)fv  (y)dy  .       (3) 
-co   R,Yp  _m  p        Yp 


The  right  side  of  equation  (3)  is  the  expectation  of  the 
conditional  density  f(r|Y   =  y) ,  and  the  left  side  is  the 
marginal  density,  in  this  case  the  mass  density,  of  r;  hence; 


P(D)  =  p_(r)  =  /   f (r  Y  =  y)fv  (y)dy  .  (4) 

K  _oo  P  * 

Since  the  conditional  density  in  equation  (4)  was  derived 
numerically,  and  no  explicit  mathematical  expression  was 
determined,  the  integral  indicated  was  performed  by  geo- 
metric integration  techniques.   Appendix  E  shows  how  the 
area  under  the  curve 


z  =  f (r|Y  =  y)fy  (y) 


43 


was  obtained.   Equation  (4)  may  be  rewritten  as 


L/2+r 
/  f(r|Y  =y)fy  (y)dy  =  2  /      f(r|Y  =y)fy  (y)dy,   (5) 

P       P  0  P 


since 


f  (r|Y   =  y)  =  0,  \/  y?[-    L/2  -  r,  L/2  +  r]  , 


and  the  integrand  is  symmetric. 

Figure  9  is  a  plot  of  the  function  to  be  integrated; 
the  values  of  p  used  are  0,  1,  1.5,  2,    3,  5,  9,  and  °°.   The 
results  of  the  geometrical  integration  are  listed  below  in 
Table  II,  and  are  plotted  in  Figure  10. 

TABLE  II 
GEOMETRIC  INTEGRATION  RESULTS 


SHIP-SUBMARINE 

SPEED  RATIO  P (D) 

0  .9980 

1  .9415 
1.5  .7895 

2  .6816 

3  .5781 
5  .5233 
9  .5008 
00  .4885 
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Figure  9.   Plots  of  I  «  P(D)x  «y)xP(Y  »y)  vs.    y 
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Cursory  examination  of  Figure  10  reveals  that  for  speed 
ratios  less  than  four,  the  assumption  that  a  fixed  crossing 
point  will  yield  more  pessimistic  results  in  terms  of 
higher  probabilities  of  detection  than  those  arising  from 
the  assumption  of  normally-distributed  cross  points  is  not 
valid.   The  length  of  the  submarine  patrol  line  was 
varied  from  120  miles  down  to  30  miles  in  order  that  the 
effect  of  changing  the  ratio  of  patrol  line  length  to  detec- 
tion range  might  be  examined.   Figure  11  shows  the  results 
for  a  patrol  line  length  of  90  miles  (three  times  the 
detection  range) .   Note  that  the  binodal  behavior  of  the 
probability  curves  which  was  present  when  a  patrol  line 
length  of  120  miles  was  used  is  absent  in  this  case,  except 
for  the  situation  where  p  =  9  ,  where  the  binodal  behavior 
is  incipient.   The  nodal  phenomenon  observed  indicates  a 
possible  area  of  interest  for  additional  research. 
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Figure  11.   Probability  of  Detection  for  Various 
Speed  Ratios  vs.  Track  Crossing  Position 
(Perpendicular  Tracks)  for  L=90  miles 
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III.   THE  RANDOM  ENCOUNTER  -  SIMULATION 
BY  MONTE  CARLO  METHODS  - 


A.   MIDPOINT  CROSSING  MODEL 

The  encounter  model  developed  by  deterministic  methods 
was  used  to  build  the  models  and  computer  programs  for  the 
stochastic  simulation.   All  assumptions  pertaining  to  the 
deterministic  model  were  included  in  the  Monte  Carlo  ap- 
proach. 

The  submarine's  initial  position  on  his  patrol  line  was 
determined  for  each  sample  taken  by  selecting  from  a  uniform 
(0,1)  pseudorandom  number  generator.   This  procedure  agreed 
with  the  assumption  that  the  submarine  was  equally  likely  to 
be  anywhere  on  his  patrol  line  at  any  time  x. 

The  submarine's  initial  up  or  down  heading  on  its  patrol 
line  was  determined  by  selecting  a  pseudorandom  number  from 
a  uniform  (0,1)  variate  generator.   If  this  number  was  less 
than  0.5  the  submarine  was  assigned  an  initial  up  heading; 
if  the  number  was  greater  than  0.5  the  submarine  was 
assigned  an  initial  down  heading.   If  the  random  variate  was 
equal  to  0.5  another  uniform  random  variable  was  selected 
and  the  above  procedure  was  repeated. 

The  midpoint  crossing  of  the  patrol  line  by  the  merchant 
was  arranged  by  positioning  the  bottom  of  the  patrol  line  in 
such  a  manner  as  to  insure  that  this  phenomenon  occurred. 
Section  (11)  gives  the  analytical  explanation  dealing  with 
the  development  of  the  deterministic  approach. 
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The  results  obtained  using  the  Monte  Carlo  technique  are 
presented  in  Appendix  A.   Graphs  depicting  the  results  can 
be  found  in  Figure  12.   It  should  be  noted  that  the  results 
obtained  utilizing  the  stochastic  simulation  agree  with 
those  obtained  from  the  deterministic  approach. 

B.   EXTENSION  TO  NORMALLY-DISTRIBUTED  CROSSING  POINTS 

The  model  and  computer  program  for  normally-distributed 
crossing  points  of  the  convoy  route  about  the  center  of  the 
submarine  patrol  line  is  an  extension  of  the  midpoint 
crossing  model. 

The  submarine's  initial  position  and  heading  were 
determined  as  for  the  midpoint  crossing   model.   The  point 
of  intersection  of  the  submarine  patrol  line  and  the  convoy 
route  was  determined  by  generating  a  normal  N  (0,400) 
pseudorandom  variate  using  the  Box-Mueller  Method  [28] .   The 
mean  of  the  normal  variate  corresponded  to  the  center  of  the 
patrol  line.   The  problem  was  then  generated  in  a  manner 
similar  to  that  of  the  midpoint  crossing  model. 

A  relatively  small  sample  size  of  10  normal  variates, 
one  for  each  1,000  uniformly-distributed  initial  starting 
positions  for  the  submarine,  was  generated  due  to  the  time 
limitations  on  the  computer  system.   This  sample  size  took 
approximately  46  minutes  to  execute.   The  computer  program 
which  appears  in  Appendix  D  has  been  written  using  a  larger 
sample  size  in  an  effort  to  aid  in  any  future  duplications 
of  the  simulation. 
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The  results  obtained  using  this  simulation  technique 
are  tabulated  in  Appendix  A.   A  plot  of  the  results  appears 
in  Figure  13  and  Figure  14.   The  results  and  analysis  of 
this  situation  are  similar  to  those  obtained  from  the  method 
of  analysis  by  numerical  integration  in  the  extension  to 
variable  crossing  points  as  is  evident  from  Figure  14. 
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IV.   AN  EXAMPLE  OF  A  NORTH  ATLANTIC  CROSSING 

A.   ASSUMPTIONS 

It  was  assumed  that  the  major  western  North  Atlantic 
seaports  are  New  York,  Norfolk,  Windward  Passage,  and 
Puerto  Rico  (Panama  Canal  traffic) ;  the  major  eastern 
North  Atlantic  seaports  were  assumed  to  be  the  Faero  Islands 
(representing  the  straits  between  Iceland  and  Scotland) , 
the  English  Channel,  and  Gibraltar.   This  assumption  results 
in  twelve  possible  Atlantic-crossing  convoy  routes.   These 
possible  convoy  routes  were  used  as  a  basis  to  determine 
the  following  assumptions: 

1)  allocation  of  enemy  submarines 

2)  lengths  of  submarine  barrier  patrol  lines 

3)  angle  of  orientation  of  the  submarine  patrol  lines 
to  the  convoy  tracks. 

A  consequence  resulting  from  the  assumption  of  specific 
convoy  routes  is  a  deviation  from  the  assumptions  made  by 
Koopman  [18],  Neufer  [5],  Dobbie  [24],  and  others  which  was 
that  of  a  uniform  distribution  of  submarines  throughout  an 
area.   Figure  15  is  a  sketch  which  illustrates  the  rationale 
for  the  assumptions  made  herein. 

Given  that  the  only  possible  convoy  routes  are  as  shown 
in  the  figure,  an  enemy  might  be  expected  to  allocate  his 
anti-shipping  submarine  patrols  as  follows: 

1.   One  submarine  to  patrol  each  port  (since  ASW 
activity  would  be  expected  to  be  heaviest  near  these  ports, 
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greater  numbers  of  patrols  at  these  locations  was  not 
assumed) . 

2.  A  submarine  patrol  at  each  expected  juncture  of 
convoy  routes  (this  effectively  places  a  barrier  on  each 
of  the  convoy  routes  which  merge  at  the  junction) . 

3.  No  segment  of  convoy  track  greater  than  50  0  miles 
to  be  unopposed  by  a  barrier. 

The  preceding  assumptions  require  that  an  enemy  have 
the  capability  to  maintain  70  submarines  on  station  in 
the  North  Atlantic.   Note  that  this  distribution  results 
in  105  barriers  against  the  possible  convoy  routes  assumed, 
even  though  only  70  submarines  are  employed. 

Another  assumption  made  was  that  the  submarine  barriers 
would  be  oriented  according  to  the  following  rules: 

1.  Patrolling  of  straits  to  be  oriented  so  as  to 
minimize  the  patrol  line  length. 

2.  Patrolling  of  seaport  exits  (entrances)  to  be  per- 
pendicular to,  and  centered  on,  a  bisector  of  the  widest 
angle  formed  by  two  departure  courses. 

3.  Patrolling  of  route  junctions  to  be  along  the 
bisector  of  the  largest  angle  formed  by  the  intersection 
of  the  two  convoy  routes ,  and  to  be  centered  on  the 
expected  junction. 

4.  All  other  patrols  to  be  centered  on,  and  oriented 
at  right  angles  to  the  expected  route. 

Additionally,  all  submarine  patrolling  speeds  were  assumed 
to  be  10  knots,  and  all  patrol  lines  were  assumed  to  be 
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120  miles  in  length  (unless  shortened  by  geographical 
limitations) . 

B.   APPLICATION  OF  THE  MODELS 

1 .   Orientation  of  Barriers  on  the  Norfolk-English 
Channel  Route 

The  route  between  Norfolk  and  the  English  Channel 

is  approximately  3,500  miles  long,  and  intersects  the  routes 

from  New  York  to  Gibraltar,  Windward  Passage  to  the  Faero 

Islands,  and  Puerto  Rico  to  the  Faero  Islands.   The  rules 

assumed  required  that  the  barrier  patrols  on  this  convoy 

route  be  given  the  following  parameters: 

TABLE  III 

PARAMETERS  OF  THE  BARRIERS  FOR  THE  EXAMPLE 

BARRIER 

#  BARRIER  LOCATION         a°        L 

1  Norfolk 

2  N.Y. -Gibraltar  junction 
3,4       Along  the  route 

5  W. P. -Faero  junction 

6  P. R. -Faero  junction 

7  English  Channel 

(Land's  End-Ushant  Is.)      0      90 


5 

120  miles 

12 

120 

0 

120 

10 

120 

15 

120 

2.   So 


lution  to  the  Example  Problem 


In  general, 

P (Detection)  =  1  -  P (No  Detection), 
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and  in  this  case, 

P (No  Detection)  =  P (No  Detection  by  Barrier  #1) 

xP (No  Detection  by  Barrier  #2) 

J\.      •  •  •   X 

xP  (No  Detection  by  Barrier  n)  , 
assuming  stochastic  independence  for  simplicity.   In  this 
example,  n  =  7.   Since  the  probability  of  detection  is  a 
function  of  the  angle  of  track  a,  the  patrol  line  length  L, 
the  detection  range  r,  and  the  ship-submarine  speed  ratio 
p,  the  probability  of  detection  for  barrier  i  may  be 
represented  as 


Pi(D)  =  P(D;ai,Li,ri,pi)   , 


and 


P(D)  =   n  [1  -  P(D;ai#Li,ri,pi)] 
i=l 


According  to  the  assumptions  made, 


rl  :  r2    *  *  *    r7  :  r 


P-j_  =  P2  =  •  •  •  =  P7  =  P 


/ 


and, 


P  (Detection)  =   II  [1  -  P  (D;  a.  ,L.  ,r ,  p)  ]  (6) 

i=l  1   1 


The  probabilities  of  detection  for  each  of  the  barriers 
listed  in  Table  3  are  tabulated  in  Table  IV. 
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Using  the  results  from  Table  IV  in  equation  (6), 

1  -  P(D|p=p*)  =  [1  -  P(D|p=p*,L  =  120, a  =  0,r  =  30)]6 

x[l  -  P(D|p=p*,L  =  90,a=0,r  =  30)] 
The  results  are  tabulated  in  Table  5. 

TABLE  V 

P (DETECTION)  VERSUS  SPEED  FOR 
THE  ROUTE  OF  THE  EXAMPLE 


p  *  P  (D| p=p*) 

0  1.0000 

1  1.0000 
1.5  1.0000 

2  1.0000 
2.5  1.0000 

3  .9996 

4  .9986 

5  .9978 
9  .9957 

.9910 


Table  V  indicates,  ceterus  paribus,  that  speed  advantage 
alone  for  the  surface  ship  will  not  lessen  the  probability 
of  detection  significantly  by  patrolling  submarines  for  a 
typical  North  Atlantic  crossing. 
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V.   SUMMARY  AND  CONCLUSIONS 

A.  MONTE  CARLO  VERSUS  NUMERICAL  INTEGRATION  METHODS 
Figure  14  shows  that  close  agreement  exists  between  the 

Monte  Carlo  and  numerical  integration  techniques;  the  numer- 
ical integration  method  is  more  economical  in  terms  of 
computer  usage  time. 

B.  THE  INDIVIDUAL  RANDOM  ENCOUNTER 

1.  Midpoint  Crossing  at  Various  Angles 

Detection  probability  decreased  from  unity  at  P  i  1 
for  all  angles  of  intersection  between  the  ship's  track  and 
the  submarine  barrier  patrol  line,  to  0.5  at  p  >.  9  (for 
intersection  at  right  angles) .   Ninety- five  per  cent  of  the 
decrease  occurred  at  speed  ratios  of  eight  and  below  for 
all  angles,  and  for  speed  ratios  of  4.5  and  below  for  inter- 
sections at  right  angles.   Stated  another  way,  for  assumed 
submarine  speed  of  10  knots,  and  using  detection  probability 
as  the  sole  criterion  for  selecting  ship  speed  capabilities, 
merchant  ship  speeds  greater  than  40  knots  do  not  seem 
warranted. 

2 .  Normally-Distributed  Crossing  Points 

a.   The  assumption  of  normally-distributed  crossing 
points  resulted  in  lower  detection  probabilities  for  ship- 
submarine  speed  ratios  less  than  four,  than  did  the  assump- 
tion of  a  fixed  crossing  point  at  the  center  of  the  patrol 
line.   For  speed  ratios  greater  than  four,  both  assumptions 
yielded  essentially  identical  results. 
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b.   If  ship-submarine  speed  ratios  equal  to  or 
greater  than  four  are  to  be  considered,  then  the  binodal 
behavior  of  the  probabilities  of  detection,  illustrated  in 
Figure  8,  should  be  examined  further.   It  appears  that 
this  behavior  may  also  be  a  function  of  the  ratio  of  sub- 
marine barrier  patrol  length  to  assumed  detection  range. 

C.   NORTH  ATLANTIC  TRANSIT  EXAMPLE 

When  six  or  more  submarine  barriers  are  crossed  by  the 
merchant  ship's  track,  the  reduction  in  P (Detection) 
resulting  from  increasing  the  merchant  ship's  speed  is 
negligible.   It  is  evident  that,  from  the  viewpoint  of 
decreasing  random  detection  probability,  a  large  speed 
advantage  is  not  the  answer.   If  a  large  speed  advantage 
by  a  merchant  ship  over  a  patrolling  submarine  increases 
the  survivability  of  that  merchant  ship,  then  the  increase 
in  survivability  is  due,  not  to  decreased  detections,  but 
to  other  factors,  e.g.,  reduced  submarine  weapon  system 
effectiveness  or  reduced  exposure  time. 
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APPENDIX  A 


PROBABILITY  CF  RANDCM  ENCOUNTER 
BASIC  MODEL 

OcFINITE  RANGE  LAW 

DFTFCTICN  RANGE:   30. C  MILES 

LENGTH  OF  SL8  PATROL  LINE:  120.0  *ILES 

X-DISTANCE  TO  PATROL  LINE:   30.0  *ILES 

SUBMARINE  SPEED:  1G.0  KNOTS 

NO.  PF  SAMPLES  TAKEN  FOR  EACH  SCENARIO:   96C 

AVC.NC.  SHIP   SPEED   RELATIVE  SPEEDS    ROUTE  PPCBABILITY  OF 
LEGS    SPEfcP  PATIO   UPLEG    DCWNLEG    ANGLE     DETECTION 


1.4C 

ICC 

1.0 

14.1 

14.1 

0 

1.0000 

1.36 

15.0 

1.5 

18.0 

18.0 

0 

0.9354 

1.3? 

2  CO 

2.0 

22.4 

22.4 

0 

C.8C00 

1.2P 

25.0 

2.5 

26.9 

26.9 

0 

0.6937 

1.25 

3CC 

3.0 

31.6 

31.6 

0 

0.6250 

1.2? 

35.0 

3.5 

36.4 

36.4 

C 

0.5750 

1.21 

40.0 

4.0 

41.2 

41.2 

0 

0.5437 

1.18 

5C.C 

5.0 

51. C 

51.0 

0 

0.5104 

1.17 

6C0 

6.0 

6C.8 

60.8 

0 

0.5062 

1.15 

7C.0 

7.0 

70.7 

70.7 

0 

0.5062 

1.15 

ec.o 

CO 

BO. 6 

80.6 

0 

0.5021 

1.14 

90.  C 

9.0 

9C.6 

90.6 

0 

0.5021 

1.39 

ICO 

1.0 

10.  7 

16.9 

25 

1.0000 

1.36 

15.0 

1.5 

14.1 

21.3 

25 

0.9427 

1.32 

2C.C 

2.0 

18.2 

25.9 

25 

C.8698 

1.29 

25.0 

2.5 

22.7 

30.6 

25 

0.7781 

1.27 

3  CO 

3.0 

27.3 

35.4 

25 

0.7240 

1.2* 

35.0 

3.5 

32.1 

40.3 

25 

0.671-5 

1.2  2 

4C.C 

*.o 

36.9 

45.1 

25 

0.6406 

L.19 

r-c.o 

5.0 

46.7 

55.0 

25 

C.5865 

1.17 

ICC 

6.0 

56.5 

64.9 

25 

0.5687 

1.16 

7C.C 

7.0 

66.4 

74.8 

25 

0.5562 

1.15 

80.0 

CO 

76.3 

84.7 

25 

0.5542 

1.1* 

9C.0 

9.0 

86.3 

94.7 

25 

0.5542 
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AVG.NO.  ship  speed 

LEGS    SPEED  RATIO 


RELATIVE  SPEEDS 
UPLEG    DOWNLEG 


ROUTE  PROBABILITY  OF 
ANGLE     DETECTION 


1.39 

10.0 

1.0 

10. C 

17.3 

30 

1.0000 

1.36 

15.0 

1.5 

13.2 

21.8 

30 

0.9469 

1.32 

20.0 

2.0 

17.3 

26.5 

30 

0.8833 

1.28 

25.0 

2.5 

21.8 

31.2 

30 

0.8187 

1.26 

30.0 

3.0 

26.5 

36.1 

30 

0.7615 

1.25 

35.0 

3.5 

31.2 

40.9 

30 

0.7229 

1.23 

40.0 

4.0 

36.1 

45.8 

30 

0.6865 

1.20 

50.0 

5.0 

45.8 

55.7 

30 

0.6344 

1.18 

60.0 

6.0 

55.7 

65.6 

30 

0.6052 

1.16 

7C.0 

7.0 

65.6 

75.5 

30 

0.5906 

1.15 

80.0 

8.0 

75.5 

85.4 

30 

0.5802 

1.15 

90.0 

9.0 

85.4 

95.4 

30 

0.5792 

1.39 

10.0 

1.0 

9.2 

17.7 

35 

1.0000 

1.35 

15.0 

1.5 

12.4 

22.3 

35 

0.9552 

1.31 

20.0 

2.0 

16.4 

27.0 

35 

0.8969 

1.28 

25.0 

2.5 

20.9 

31.8 

35 

0.8583 

1.26 

30.0 

3.0 

25.6 

36.7 

35 

0.8083 

1.24 

35.0 

3.5 

30.4 

41.6 

35 

0.7677 

1.23 

40.0 

4.0 

35.2 

46.5 

35 

0.7396 

1.20 

50.0 

5.0 

45.0 

56.3 

35 

0.6948 

1.18 

60.0 

6.0 

54.9 

66.2 

35 

0.6583 

1.17 

70.0 

7.0 

64.8 

76.2 

35 

0.6375 

1.16 

80.0 

8.0 

74.7 

86.1 

35 

0.6250 

1.15 

90.0 

9.0 

84.7 

96.1 

35 

0.6177 

1.39 

10.0 

1.0 

8.5 

18.1 

40 

1.0000 

1.34 

15.0 

1.5 

11.5 

22.8 

40 

0.9677 

1.31 

20.0 

2.0 

15.6 

27.5 

40 

0.9177 

1.28 

25.0 

2.5 

20.1 

32.3 

40 

0.8812 

1.25 

30.0 

3.0 

24.8 

37.2 

40 

0.8594 

1.23 

35.0 

3.5 

29.6 

42.1 

40 

0.8323 

1.22 

40.0 

4.0 

34.4 

47.1 

40 

0.8021 

1.20 

50.0 

5.0 

44.2 

56.9 

40 

0.7573 

1.19 

60.0 

6.0 

54.1 

66.9 

40 

0.7323 

1.17 

70.0 

7.0 

64.0 

76.8 

40 

0.7052 

1.16 

80.0 

8.0 

74.0 

86.8 

40 

0.6906 

1.15 

90.0 

9.0 

83.9 

96.7 

40 

0.6750 
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VG.NO. 

SHIP 

SPEED 

RELATIVE 

SPEEDS 

ROUTE 

PROBABILITY 

LEGS 

SPEtD 

RATIO 

UPLEG 

DOWNLEG 

ANGLE 

DETECTIOI 

1.39 

10.0 

1.0 

7.7 

18.5 

45 

1.0000 

1.3* 

15.0 

1.5 

10.6 

23.2 

45 

0.9865 

1.3C 

20.0 

2.0 

14.7 

28.0 

45 

0.9365 

1.27 

25.0 

2.5 

19.3 

32.9 

45 

0.9115 

1.25 

30.0 

3.0 

24.0 

37.7 

45 

0.8927 

1.23 

35.0 

3.5 

28.8 

42.7 

45 

0.8802 

1.21 

4C.C 

4.0 

33.7 

47.6 

45 

0.87C8 

1.19 

50.0 

5.0 

43.5 

57.5 

45 

0.8344 

1.18 

6C.0 

6.0 

53.4 

67.4 

45 

0.8073 

1.17 

70.0 

7.0 

63.3 

77.4 

45 

0.7875 

1.16 

80. C 

8.0 

73.3 

87.4 

45 

0.7740 

1.16 

°0.0 

9.0 

83.2 

97.3 

45 

0.7625 

1.39 

10.0 

1.0 

6.8 

18.8 

50 

1.0000 

1.33 

15.0 

1.5 

9.8 

23.6 

50 

l.OOCO 

1.28 

20.0 

2.0 

13.9 

28.4 

50 

0.9656 

1.25 

25.0 

2.5 

18.5 

33.3 

50 

0.9427 

1.23 

3C.C 

3.0 

23.2 

38.2 

50 

0.9302 

1.21 

35.0 

3.5 

28.1 

43.1 

50 

0.9240 

1.2C 

40. 0 

4.0 

33.0 

48.1 

50 

0.9177 

1.19 

50.0 

5.0 

42.8 

58.0 

50 

0.9083 

1.17 

60.0 

6.0 

52.7 

68.0 

50 

0.9021 

1.16 

7C.0 

7.0 

62.7 

77.9 

50 

0.8854 

1.15 

8C.0 

8.0 

72.6 

87.9 

50 

0.8698 

1.14 

90. C 

9.0 

82.6 

97.9 

50 

0.8583 

1.39 

10. c 

1.0 

6.C 

19.1 

55 

1.0000 

1.31 

15.0 

1.5 

8.9 

23.9 

55 

1.0000 

1.26 

20. 0 

2.0 

13.1 

28.8 

55 

1.0000 

1.22 

25.0 

2.5 

17.8 

33.7 

55 

0.9844 

1.20 

30.0 

3.0 

22.6 

38.6 

55 

0.9719 

1.16 

35.0 

3.5 

27.4 

43.6 

55 

0.9656 

1.17 

40.  C 

4.0 

32.3 

48.5 

55 

0.9615 

1.15 

5C.0 

5.0 

42.2 

58.5 

55 

0.9552 

1.13 

60.0 

6.0 

52.1 

68.4 

55 

0.9510 

1.12 

70. C 

7.0 

62.1 

78.4 

55 

0.9490 

1.12 

80. G 

8.0 

72.0 

88.4 

55 

0.9469 

1.11 

90.0 

°.0 

82.0 

98.4 

55 

0.9448 
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■     :■ 

OF 


AVG.NO. 

SHIP 

SPEED 

RELATIVE 

SPEEDS 

ROUTE 

PROBABILITY 

LEGS 

SPEED 

RATIO 

UPLEG 

OOWNLEG 

ANGLE 

DETECTIO 

1.39 

ICO 

1.0 

5.2 

19.3 

60 

1.0000 

1.29 

15.0 

1.5 

8.1 

24.2 

60 

1.0000 

1.23 

20.0 

2.0 

12.4 

29.1 

60 

1.0000 

1.19 

25.0 

2.5 

17.1 

34.0 

60 

1.0000 

1.16 

30.0 

3.0 

21.9 

39.0 

60 

1.0000 

1.14 

35.0 

3.5 

26.8 

43.9 

60 

.1.0000 

1.12 

4C.0 

4.0 

31.7 

48.9 

60 

1.0000 

I.  10 

50.0 

5.0 

41.6 

58.9 

60 

1.0000 

1.C8 

60.0 

6.0 

51.6 

68.8 

60     ' 

1.0000 

I.C7 

70.  C 

7.0 

61.5 

78.8 

60 

1.0000 

1.C6 

80.0 

e.o 

71.5 

88.8 

60 

1.0000 

1.C5 

90.0 

9.0 

81.5 

98.8 

60 

1.0000 
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PROBABILITY  CF  RANDOM  ENCOUNTER 
EXTENDED  MODEL 


DEFINITF  RANGE  LAW 

DETECTION  RANGE:   30.0  MILES 

LENGTH  OF  SUB  PATRCL  LINE:  120.0  MILES 

X-DISTANCE  TO  PATRCL  LINE:   30.0  MILES 

CROSSING  ANGLE:  CO  CEGREES  (PERPENDICULAR) 

SUBMARINE  SPEED:  10.0  KNCTS 

NO.  CF  SAMPLES  TAKEN  FCR  EACH  SCENARIO:   960 


AVG.NO.  SHIP   SPEED 
LEGS    SPEED  RATIO 


RELATIVE  SPEEDS 
UPLEG    DOWNLEG 


ROUTE  PROBABILITY  OF 
POSIT     DETECTION 


1.40 

ICO 

1.0 

14.1 

14.1 

60.0 

l.OOOC 

1.36 

15. C 

1.5 

18.0 

18.0 

60.0 

0.9354 

1.33 

2C.C 

2.0 

22.4 

22.4 

60.0 

ceooo 

1.28 

25.0 

2.5 

26.9 

26.9 

60.0 

0.6937 

1.25 

3CC 

3.0 

31.6 

31.6 

60.0 

0.6250 

1.21 

40. C 

4.0 

41.2 

41.2 

60.0 

0.5437 

1.18 

50.0 

5.0 

51.0 

51.0 

60.0 

0.5104 

2.32 

9C.6 

9.0 

90.6 

90.6 

60.0 

0.5021 

1.4C 

ICG 

1.0 

14.1 

14.1 

65.0 

1.0000 

1.36 

15.0 

1.5 

18.0 

18.0 

65.0 

0.926C 

1.32 

2C.0 

2.0 

22.4 

22.4 

65.0 

0.7844 

1.28 

25.0 

2.5 

26.9 

26.9 

65.0 

0.6937 

1.25 

30. C 

3.0 

31.6 

31.6 

65.0 

0.6250 

1.21 

*CC 

4.0 

41.2 

41.2 

65.0 

0.5500 

1.18 

5C.0 

5.0 

51.0 

51.0 

65.3 

0.5260 

2.32 

90.0 

9.0 

90.6 

90.6 

65.0 

0.5021 

1.40 

ICO 

1.0 

14.1 

14.1 

70.0 

l.OOOC 

1.35 

15.0 

1.5 

18.0 

18.0 

7C.0 

0.8687 

1.30 

2C0 

2.0 

22.4 

22.4 

70.0 

0.7635 

1.26 

25.0 

2.5 

26.9 

26.9 

70.0 

0.6750 

1.24 

30.0 

3.0 

31.6 

31.6 

70.0 

0.6177 

1.21 

40.0 

4.0 

41.2 

41.2 

70.0 

0.5708 

1.18 
2.32 

50.0 

5.0 

51.0 

51.0 

70.0 

0.5469 

9C.C 

9.0 

90.6 

90.6 

70.0 

0.5021 

1.40 

10.0 

1.0 

14.1 
16.0 

14.1 

75.0 

0.9781 

1.33 

15.0 

1.5 

18.0 

75.0 

0.8062 

1.28 

2C.0 

2.0 

22.4 

22.4 

75.0 

0.7125 

1.2* 

25. C 

2.5 

26.9 

26.9 

75.0 

0.6542 

1.22 

3C.0 

3.0 

31.6 

31.6 

75.0 

0.6104 

1.19 

40.0 

4.0 

41.2 

41.2 

75.0 

0.5771 

i.18 

5C.0 

5.C 

51.0 

51.0 

75.0 

0.5604 

2.31 

9C.0 

9.0 

90.6 

90.6 

75.0 

0.5229 
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AVG.NO.  ship  speed 
LEGS    SPEED  RATIO 


RELATIVE  SPEEDS 
UPLEG    DOWNLEG 


RCUTE  PROBABILITY  OF 
POSIT     DETECTION 


1.40 

10.0 

1.0 

14.1 

14.1 

80.0 

0.9365 

1.31 

15.0 

1.5 

18.0 

18.0 

80.0 

0.7437 

1.26 

20.0 

2.0 

22.4 

22.4 

80.0 

0.6500 

1.22 

25.0 

2.5 

26.9 

26.9 

80.0 

0.6031 

1.20 

30.0 

3.0 

31.6 

31.6 

80.0 

0.5969 

1.17 

40.0 

4.0 

41.2 

41.2 

80.0 

0.5771 

1.16 

50. C 

5.0 

51.0 

51.0 

80.0 

0.5604 

2.28 

90.0 

9.0 

90.6 

90.6 

80.0 

0.5292 

1.39 

10.0 

1.0 

14.1 

14.1 

85.0 

0.8844 

1.29 

15.0 

1.5 

18.0 

18.0 

85.0 

0.6812 

1.23 

20.0 

2.0 

22.4 

22.4 

85.0 

0.5875 

1.20 

25.0 

2.5 

26.9 

26.9 

85.0 

0.5615 

1.18 

30.0 

3.0 

31.6 

31.6 

85.0 

0.5552 

1.15 

40.0 

4.0 

41.2 

41.2 

85.0 

0.5490 

1.13 

50.0 

5.0 

51.0 

51.0 

85.0 

0.5469 

2.24 

90.0 

9.0 

90.6 

90.6 

85.0 

0.5292 

1.36 

10.0 

1.0 

14.1 

14.1 

90.0 

0.8219 

1.27 

15.0 

1.5 

18.0 

18.0 

90.0 

0.6187 

1.21 

20.0 

2.0 

22.4 

22.4 

90.0 

0.5302 

1.18 

25.0 

2.5 

26.9 

26.9 

90.0 

0.5198 

1.16 

30.0 

3.0 

31.6 

31.6 

90.0 

0.5135 

1.13 

40.0 

4.0 

41.2 

41.2 

90.0 

0.5073 

1.11 

50.0 

5.0 

51.0 

51.0 

90.0 

0.5052 

2.19 

90.0 

9.0 

90.6 

90.6 

90.0 

0.5010 

1.30 

10.0 

1.0 

14.1 

14.1 

110.0 

0.5719 

1.21 

15.0 

1.5 

18.0 

18.0 

110.0 

0.3844 

1.17 

20.0 

2.0 

22.4 

22.4 

110.0 

0.3635 

1.14 

25.0 

2.5 

26.9 

26.9 

110.0 

0.3531 

1.13 

30.0 

3.0 

31.6 

31.6 

110.0 

0.3469 

1.11 

40.0 

4.0 

41.2 

41.2 

110.0 

0.3406 

1.09 

50.0 

5.0 

51.0 

51.0 

110.0 

0.3385 

2.16 

90.0 

9.0 

90.6 

90.6 

110.0 

0.3344 

1.33 

10.0 

1.0 

14.1 

14.1 

130.0 

0.3219 

1.23 

15.0 

1.5 

18.0 

18.0 

130.0 

0.2177 

1.18 

20.0 

2.0 

22.4 

22.4 

130.0 

0.1969 

1.14 

25.0 

2.5 

26.9 

26.9 

130.0 

0.1865 

1.13 

30. 0 

3.0 

31.6 

31.6 

130.0 

0.1802 

1.11 

40.0 

4.0 

41.2 

41.2 

130.0 

0.1740 

1.09 

50.0 

5.0 

51.0 

51.0 

130.0 

0.1719 

2.16 

90.0 

9.0 

90.6 

90.6 

130.0 

0.1677 

1.58 

ICO 

1.0 

14.1 

14.1 

150.0 

0.0 

1.42 

15.0 

1.5 

18.0 

18.0 

150.0 

0.0 

1.33 

20.0 

2.0 

22.4 

22.4 

150.0 

0.0010 

1.28 

25.0 

2.5 

26.9 

26.9 

150.0 

0.0010 

1.25 

30.0 

3.0 

31.6 

31.6 

150.0 

0.0010 

1.21 

40.0 

4.0 

41.2 

41.2 

150.0 

0.0010 

1.18 

50.0 

5.0 

51.0 

51.0 

150.0 

0.0010 

2.32 

9C.0 

9.0 

90.6 

90.6 

150.0 

0.0 
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MONTE  CARLO  TECHNIQUE 
PROBABILITY  OF  RANDOM  ENCOUNTER 
MID-POINT  CROSSING 


DEFINITE  RANGE  LAW 

DETECTION  RANGE:   30.0  MILES 

LENGTH  OF  SUB  PATROL  LINE:   120.0  MILES 

X-DISTANCE  TO  PATROL  LINE:   30.0  MILES 

SUBMARINE  SPEED:   10.0  KNOTS 

NO.  OF  SAMPLES  TAKEN  FOR  EACH  SCENARIO:   1000 


AVG.  NO. 

SHIP 

SPEED 

RELATIVE  SPEED 

ALFA 

PROB. 

LEGS 

SPEED 

RATIO 

UPLEG 

DOWNLEG 

DETECTION 

1.41 

10.0 

1.0 

14.1 

14.1 

0.0 

1.0000 

1.34 

20.0 

2.0 

22.4 

22.4 

0.0 

0.8200 

1.25 

30.0 

3.0 

31.6 

31.6 

0.0 

0.6270 

1.21 

40.0 

4.0 

41.2 

41.2 

0.0 

0.5470 

1.18 

50.0 

5.0 

51.0 

51.0 

0.0 

0.5290 

1.17 

60.0 

6.0 

60.8 

60.8 

0.0 

0.5060 

1.15 

70.0 

7.0 

70.7 

70.7 

0.0 

0.4810 

1.16 

80.0 

8.0 

80.6 

80.6 

0.0 

0.4810 

1.13 

90.0 

9.0 

90.6 

90.6 

0.0 

0.5070 

1.42 

10.0 

1.0 

13.5 

14.7 

0.0873 

1.0000 

1.32 

20.0 

2.0 

21.6 

23.1 

0.0873 

0.7650 

1.25 

30.0 

3.0 

30.8 

32.4 

0.0873 

0.6420 

1.19 

40.0 

4.0 

40.4 

42.1 

0.0873 

0.5460 

1.18 

50.0 

5.0 

50.1 

51.8 

0.0873 

0.5210 

1.19 

60.0 

6.0 

60.0 

61.7 

0.0873 

0.5020 

1.15 

70.0 

7.0 

69.8 

71.6 

0.0873 

0.4980 

1.14 

80.0 

8.0 

79.8 

81.5 

0.0873 

0.5180 

1.14 

90.0 

9.0 

89.7 

91.4 

0.0873 

0.4930 

1.38 

10.0 

1.0 

12.9 

15.3 

0.1745 

1.0000 

1.34 

20.0 

2.0 

20.7 

23.9 

0.1745 

0.8140 

1.24 

30.0 

3.0 

29.9 

33.2 

0.1745 

0.6380 

1.23 

40.0 

4.0 

39.5 

42.9 

0.1745 

0.5640 

1.17 

50.0 

5.0 

49.3 

52.7 

0.1745 

0.5190 

1.15 

60.0 

6.0 

59.1 

62.5 

0.1745 

0.5130 

1.15 

70.0 

7.0 

69.0 

72.4 

0.1745 

0.5210 

1.14 

80.0 

8.0 

78.9 

82.3 

0.1745 

0.5070 

1.14 

90.0 

9.0 

88.8 

92.3 

0.1745 

0.5000 
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AVG.  NO. 

SHIP 

SPEED 

RELATIVE  SPEED 

ALFA 

PROB. 

LEGS 

SPEED 

RATIO 

UPLEG 

DOWNLEG 

DETECTION 

1.40 

10.0 

1.0 

12.2 

15.9 

0.2618 

1.0000 

1.32 

20.0 

2.0 

19.9 

24.6 

0.2618 

0.8120 

1.26 

30.0 

3.0 

29.1 

34.0 

0.2618 

0.6500 

1.22 

40.0 

4.0 

38.6 

43.7 

0.2618 

0.5530 

1.20 

50.0 

5.0 

48.4 

53.5 

0.2618 

0.5260 

1.17 

60.0 

6.0 

58.2 

63.3 

0.2618 

0.5270 

1.15 

70.0 

7.0 

68.1 

73.2 

0.2618 

0.5190 

1.14 

80.0 

8.0 

78.0 

83.2 

0.2618 

0.5150 

1.15 

90.0 

9.0 

87.9 

93.1 

0.2618 

0.5190 

1.39 

10.0 

1.0 

11.5 

16.4 

0.3491 

1.0000 

1.32 

20.0 

2.0 

19.1 

25.2 

0.3491 

0.8600 

1.24 

30.0 

3.0 

28.2 

34.7 

0.3491 

0.6920 

1.20 

40.0 

4.0 

37.8 

44.4 

0.3491 

0.6070 

1.19 

50.0 

5.0 

47.5 

54.2 

0.3491 

0.5470 

1.18 

60.0 

6.0 

57.4 

64.1 

0.3491 

0.5480 

1.16 

70.0 

7.0 

67.2 

74.0 

0.3491 

0.5290 

1.14 

80.0 

8.0 

77.2 

83.9 

0.3491 

0.5340 

1.16 

90.0 

9.0 

87.1 

93.9 

0.3491 

0.5070 

1.42 

10.0 

1.0 

10.7 

16.9 

0.4363 

1.0000 

1.30 

20.0 

2.0 

18.2 

25.9 

0.4363 

0.8570 

1.26 

30.0 

3.0 

27.3 

35.4 

0.4363 

0.7380 

1.22 

40.0 

4.0 

36.9 

45.1 

0.4363 

0.6430 

1.18 

50.0 

5.0 

46.7 

55.0 

0.4363 

0.5910 

1.18 

60.0 

6.0 

56.5 

64.9 

0.4363 

0.5530 

1.15 

70.0 

7.0 

66.4 

74.8 

0.4363 

0.5860 

1.16 

80.0 

8.0 

76.3 

84.7 

0.4363 

0.5360 

1.13 

90.0 

9.0 

86.3 

94.7 

0.4363 

0.5490 

1.41 

10.0 

1.0 

10.0 

17.3 

0.5236 

1.0000 

1.30 

20.0 

2.0 

17.3 

26.5 

0.5236 

0.8890 

1.26 

30.0 

3.0 

26.5 

36.1 

0.5236 

0.7270 

1.24 

40.0 

4.0 

36.1 

45.8 

0.5236 

0.6690 

1.21 

50.0 

5.0 

45.8 

55.7 

0.5236 

0.6100 

1.18 

60.0 

6.0 

55.7 

65.6 

0.5236 

0.5830 

1.18 

70.0 

7.0 

65.6 

75.5 

0.5236 

0.5910 

1.15 

80.0 

8.0 

75.5 

85.4 

0.5236 

0.5780 

1.13 

90.0 

9.0 

85.4 

95.4 

0.5236 

0.5540 

1.37 

10.0 

1.0 

9.2 

17.7 

0.6109 

1.0000 

1.31 

20.0 

2.0 

16.4 

27.0 

0.6109 

0.9010 

1.26 

30.0 

3.0 

25.6 

36.7 

0.6109 

0.8140 

1.25 

40.0 

4.0 

35.2 

46.5 

0.6109 

0.7480 

1.19 

50.0 

5.0 

45.0 

56.3 

0.6109 

0.7090 

1.19 

60.0 

6.0 

54.9 

66.2 

0.6109 

0.6540 

1.17 

70.0 

7.0 

64.8 

76.2 

0.6109 

0.6190 

1.16 

80.0 

8.0 

74.7 

86.1 

0.6109 

0.6290 

1.14 

90.0 

9.0 

84.7 

96.1 

0.6109 

0.6310 
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G.  NO. 

SHIP 

SPEED 

RELATI 

VE  SPEED 

ALFA 

PROB. 

LEGS 

SPEED 

RATIO 

UPLEG 

DOWNLEG 

DETECTION 

1.40 

10.0 

1.0 

8.5 

18.1 

0.6981 

1.0000 

1.31 

20.0 

2.0 

15.6 

27.5 

0.6981 

0.9350 

1.24 

30.0 

3.0 

24.8 

37.2 

0.6981 

0.8720 

1.23 

40.0 

4.0 

34.4 

47.1 

0.6981 

0.8090 

1.22 

50.0 

5.0 

44.2 

56.9 

0.6981 

0.7590 

1.19 

60.0 

6.0 

54.1 

66.9 

0.6981 

0.7560 

1.16 

70.0 

7.0 

64.0 

76.8 

0.6981 

0.7080 

1.16 

80.0 

8.0 

74.0 

86.8 

0.6981 

0.7080 

1.14 

90.0 

9.0 

83.9 

96.7 

0.6981 

0.6880 

1.42 

10.0 

1.0 

7.7 

18.5 

0.7854 

1.0000 

1.30 

20.0 

2.0 

14.7 

28.0 

0.7854 

0.9410 

1.25 

30.0 

3.0 

24.0 

37.7 

0.7854 

0.8740 

1.20 

40.0 

4.0 

33.7 

47.6 

0.7854 

0.8550 

1.20 

50.0 

5.0 

43.5 

57.5 

0.7854 

0.8310 

1.20 

60.0 

6.0 

53.4 

67.4 

0.7854 

0.7850 

1.16 

70.0 

7.0 

63.3 

77.4 

0.7854 

0.8060 

1.16 

80.0 

8.0 

73.3 

87.4 

0.7854 

0.7680 

1.17 

90.0 

9.0 

83.2 

97.3 

0.7854 

0.7560 

1.38 

10.0 

1.0 

6.8 

18.8 

0.8727 

1.0000 

1.26 

20.0 

2.0 

13.9 

28.4 

0.8727 

0.9570 

1.24 

30.0 

3.0 

23.2 

38.2 

0.8727 

0.9260 

1.22 

40.0 

4.0 

33.0 

48.1 

0.8727 

0.9130 

1.17 

50.0 

5.0 

42.8 

58.0 

0.8727 

0.9140 

1.17 

60.0 

6.0 

52.7 

68.0 

0.8727 

0.8940 

1.16 

70.0 

7.0 

62.7 

77.9 

0.8727 

0.8830 

1.14 

80.0 

8.0 

72.6 

87.9 

0.8727 

0.8650 

1.14 

90.0 

9.0 

82.6 

97.9 

0.8727 

0.8410 

1.39 

10.0 

1.0 

6.0 

19.1 

0.9599 

1.0000 

1.25 

20.0 

2.0 

13.1 

28.8 

0.9599 

1.0000 

1.18 

30.0 

3.0 

22.6 

38.6 

0.9599 

0.9740 

1.16 

40.0 

4.0 

32.3 

48.5 

0.9599 

0.9620 

1.16 

50.0 

5.0 

42.2 

58.5 

0.9599 

0.9520 

1.12 

60.0 

6.0 

52.1 

68.4 

0.9599 

0.9470 

1.12 

70.0 

7.0 

62.1 

78.4 

0.9599 

0.9420 

1.11 

80.0 

8.0 

72.0 

88.4 

0.9599 

0.9430 

1.11 

90.0 

9.0 

82.0 

98.4 

0.9599 

0.9410 

1.38 

10.0 

1.0 

5.2 

19.3 

1.0472 

1.0000 

1.23 

20.0 

2.0 

12.4 

29.1 

1.0472 

1.0000 

1.16 

30.0 

3.0 

21.9 

39.0 

1.0472 

1.0000 

1.12 

40.0 

4.0 

31.7 

48.9 

1.0472 

1.0000 

1.08 

50.0 

5.0 

41.6 

58.9 

1.0472 

1.0000 

1.10 

60.0 

6.0 

51.6 

68.8 

1.0472 

1.0000 

1.08 

70.0 

7.0 

61.5 

78.7 

1.0472 

1.0000 

1.07 

80.0 

8.0 

71.5 

88.8 

1.0472 

1.0000 

1.05 

90.0 

9.0 

81.5 

98.8 

1.0472 

1.0000 
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AVG.  NO. 

SHIP 

SPEED 

RELATIVE  SPEED 

ALFA 

PROB. 

LEGS 

SPEED 

RATIO 

UPLEG 

DOWNLEG 

DETECTION 

1.44 

10.0 

1.0 

4.3 

19.5 

1.1345 

1.0000 

1.24 

20.0 

2.0 

11.7 

29.4 

1.1345 

1.0000 

1.15 

30.0 

3.0 

21.4 

39.3 

1.1345 

1.0000 

1.12 

40.0 

4.0 

31.2 

49.2 

1.1345 

1.0000 

1.09 

50.0 

5.0 

41.2 

59.2 

1.1345 

1.0000 

1.08 

60.0 

6.0 

51.1 

69.2 

1.1345 

1.0000 

1.07 

70.0 

7.0 

61.1 

79.2 

1.1345 

1.0000 

1.05 

80.0 

8.0 

71.1 

89.2 

1.1345 

1.0000 

1.04 

90.0 

9.0 

81.0 

99.2 

1.1345 

1.0000 

1.55 

10.0 

1.0 

3.5 

19.7 

1.2217 

1.0000 

1.27 

20.0 

2.0 

11.1 

29.6 

1.2217 

1.0000 

1.18 

30.0 

3.0 

20.9 

39.5 

1.2217 

1.0000 

1.15 

40.0 

4.0 

30.8 

49.5 

1.2217 

1.0000 

1.11 

50.0 

5.0 

40.7 

59.5 

1.2217 

1.0000 

1.11 

60.0 

6.0 

50.7 

69.5 

1.2217 

1.0000 

1.08 

70.0 

7.0 

60.7 

79.5 

1.2217 

1.0000 

1.07 

80.0 

8.0 

70.7 

89.5 

1.2217 

1.0000 

1.07 

90.0 

9.0 

80.7 

99.5 

1.2217 

1.0000 

1.76 

10.0 

1.0 

2.6 

19.8 

1.3090 

1.0000 

1.36 

20.0 

2.0 

10.7 

29.8 

1.3090 

1.0000 

1.25 

30.0 

3.0 

20.5 

39.7 

1.3090 

1.0000 

1.18 

40.0 

4.0 

30.5 

49.7 

1.3090 

1.0000 

1.14 

50.0 

5.0 

40.4 

59.7 

1.3090 

1.0000 

1.12 

60.0 

6.0 

50.4 

69.7 

1.3090 

1.0000 

1.12 

70.0 

7.0 

60.4 

79.7 

1.3090 

1.0000 

1.08 

80.0 

8.0 

70.4 

89.7 

1.3090 

1.0000 

1.08 

90.0 

9.0 

80.4 

99.7 

1.3090 

1.0000 

2.19 

10.0 

1.0 

1.7 

19.9 

1.3963 

1.0000 

1.63 

20.0 

2.0 

10.3 

29.9 

1.3963 

1.0000 

1.44 

30.0 

3.0 

20.2 

39.9 

1.3963 

1.0000 

1.28 

40.0 

4.0 

30.2 

49.9 

1.3963 

1.0000 

1.27 

50.0 

5.0 

40.2 

59.9 

1.3963 

1.0000 

1.18 

60.0 

6.0 

50.2 

69.9 

1.3963 

1.0000 

1.15 

70.0 

7.0 

60.2 

79.9 

1.3963 

1.0000 

1.15 

80.0 

8.0 

70.2 

89.9 

1.3693 

1.0000 

1.15 

90.0 

9.0 

80.2 

99.9 

1.3963 

1.0000 
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MONTE  CARLO  TECHNIQUE 

PROBABILITY  OF  RANDOM  ENCOUNTER 

NORMALIZED  CONVOY  ROUTE 


DEFINITE  RANGE  LAW 

DETECTION  RANGE:   30.0  MILES 

LENGTH  OF  SUB  PATROL  LINE:   120.0  MILES 

X-DISTANCE  TO  PATROL  LINE:   30.0  MILES 

SUBMARINE  SPEED:   10.0  KNOTS 

NO.  OF  SAMPLES  TAKEN  FOR  EACH  SCENARIO:   1000 


AVG.  NO. 

SHIP 

SPEED 

RELATIVE  SPEED 

ALFA 

PROB. 

LEGS 

SPEED 

RATIO 

UPLEG 

DOWNLEG 

DETECTION 

1.39 

10.0 

1.0 

14.1 

14.1 

0.0 

0.9212 

1.27 

20.0 

2.0 

22.4 

22.4 

0.0 

0.6714 

1.20 

30.0 

3.0 

31.6 

31.6 

0.0 

0.5720 

1.19 

40.0 

4.0 

41.2 

41.2 

0.0 

0.5512 

1.16 

50.0 

5.0 

51.0 

51.0 

0.0 

0.5324 

1.16 

60.0 

6.0 

60.8 

60.8 

0.0 

0.5221 

1.12 

70.0 

7.0 

70.7 

70.7 

0.0 

0.5000 

1.13 

80.0 

8.0 

80.6 

80.6 

0.0 

0.5181 

1.12 

90.0 

9.0 

90.6 

90.6 

0.0 

0.5133 
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PROBABILITY  OF  RANDOM  ENCOUNTER 
EXTENDED  MODEL 


DEFINITE  RANGE  LAW 

DETECTION  RANGE:   30.0  MILES 

LENGTH  OF  SUB  PATROL  LINE:   9  0.0  MILES 

X-DISTANCE  TO  PATROL  LINE:   30.0  MILES 

CROSSING  ANGLE:   0.0  DEGREES  (PERPENDICULAR) 

SUBMARINE  SPEED:   10.0  KNOTS 

NO.  OF  SAMPLES  TAKEN  FOR  EACH  SCENARIO:   720 


AVG.  NO. 

SHIP 

SPEED 

RELATIVE  SPEED 

CROSS 

PROBABILITY 

LEGS 

SPEED 

RATIO 

UPLEG 

DOWNLEG 

POSIT 

OF  DETECTION 

1.36 

10.0 

1.0 

14.1 

14.1 

45.0 

1.0000 

1.32 

15.0 

1.5 

18.0 

18.0 

45.0 

1.0000 

1.29 

20.0 

2.0 

22.4 

22.4 

45.0 

1.0000 

1.27 

25.0 

2.5 

26.9 

26.9 

45.0 

0.9861 

1.26 

30.0 

3.0 

31.6 

31.6 

45.0 

0.9250 

1.24 

40.0 

4.0 

41.2 

41.2 

45.0 

0.8528 

1.22 

50.0 

5.0 

51.0 

51.0 

45.0 

0.8139 

2.41 

90.0 

9.0 

90.6 

90.6 

45.0 

0.7250 

1.36 

10.0 

1.0 

14.1 

14.1 

50.0 

1.0000 

1.32 

15.0 

1.5 

18.0 

18.0 

50.0 

1.0000 

1.29 

20.0 

2.0 

22.4 

22.4 

50.0 

0.9847 

1.27 

25.0 

2.5 

26.9 

26.9 

50.0 

0.9375 

1.26 

30.0 

3.0 

31.6 

31.6 

50.0 

0.9069 

1.23 

40.0 

4.0 

41.2 

41.2 

50.0 

0.8444 

1.21 

50.0 

5.0 

51.0 

51.0 

50.0 

0.7958 

2.37 

90.0 

9.0 

90.6 

90.6 

50.0 

0.7056 

1.36 

10.0 

1.0 

14.1 

14.1 

55.0 

1.0000 

1.32 

15.0 

1.5 

18.0 

18.0 

55.0 

1.0000 

1.29 

20.0 

2.0 

22.4 

22.4 

55.0 

0.9292 

1.27 

25.0 

2.5 

26.9 

26.9 

55.0 

0.8792 

1.24 

30.0 

3.0 

31.6 

31.6 

55.0 

0.8333 

1.20 

40.0 

4.0 

41.2 

41.2 

55.0 

0.7792 

1.18 

50.0 

5.0 

51.0 

51.0 

55.0 

0.7500 

2.32 

90.0 

9.0 

90.6 

90.6 

55.0 

0.7056 
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AVG .  NO . 

SHIP 

SPEED 

RELATIVE  SPEED 

CROSS 

PROBABILITY 

LEGS 

SPEED 

RATIO 

UPLEG 

DOWNLEG 

POSIT 

OF  DETECTION 

1.36 

10.0 

1.0 

14.1 

14.1 

60.0 

1.0000 

1.  32 

15.0 

1.5 

18.0 

18.0 

60.0 

0.9569 

1.28 

20.0 

2.0 

22.4 

22.4 

60.0 

0.8667 

1.24 

25.0 

2.5 

26.9 

26.9 

60.0 

0.7958 

1.21 

30.0 

3.0 

31.6 

31.6 

60.0 

0.7500 

1.17 

40.0 

4.0 

41.2 

41.2 

60.0 

0.6958 

1.15 

50.0 

5.0 

51.0 

51.0 

60.0 

0.6736 

2.26 

90.0 

9.0 

90.6 

90.6 

60.0 

0.6681 

1.36 

10.0 

1.0 

14.1 

14.1 

65.0 

1.0000 

1.32 

15.0 

1.5 

18.0 

18.0 

65.0 

0.9014 

1.26 

20.0 

2.0 

22.4 

22.4 

65.0 

0.7833 

1.22 

25.0 

2.5 

26.9 

26.9 

65.0 

0.7125 

1.19 

30.0 

3.0 

31.6 

31.6 

65.0 

0.6667 

1.16 

40.0 

4.0 

41.2 

41.2 

65.0 

0.6208 

1.14 

50.0 

5.0 

51.0 

51.0 

65.0 

0.6181 

2.23 

90.0 

9.0 

90.6 

90.6 

65.0 

0.6125 

1.36 

10.0 

1.0 

14.1 

14.1 

70.0 

1.0000 

1.31 

15.0 

1.5 

18.0 

18.0 

70.0 

0.8250 

1.24 

20.0 

2.0 

22.4 

22.4 

70.0 

0.7000 

1.21 

25.0 

2.5 

26.9 

26.9 

70.0 

0.6292 

1.18 

30.0 

3.0 

31.6 

31.6 

70.0 

0.5833 

1.15 

40.0 

4.0 

41.2 

41.2 

70.0 

0.5653 

1.13 

50.0 

5.0 

51.0 

51.0 

70.0 

0.5625 

2.23 

90.0 

9.0 

90.6 

90.6 

70.0 

0.5569 
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APPENDIX  B 
MAIN  COMPUTER  PROGRAM  FLOW  DIAGRAM 


C 


START 


) 


$L 


INITIALIZE  VARIABLES 
AND  COUNTERS 


DETERMINE  COORDINATES 

OF  LOWER  POINT  OF  SUB 

PATROL  LINE 


1801 


COMPUTE  MAX  TIME  THAT 
SHIP  WILL  BE  WITHIN 
RADIUS  OF  DETECTION 

RANGE  FROM  SUB  PATROL 


0 


/determine  initial\ 
(sub  position  vector) 
\         from  ship         / 
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10 


COMPUTE  TIME  TO  GO  FOR 
SUB  TO  REACH  END  OF 
PATROL  FROM  ITS 
INITIAL  POSITION 


COMPUTE  RELATIVE  DISTANCE 
SUB  MOVES  ON  THIS  COURSE 
(LEG) 


H 


L 


COMPUTE  FINAL  SUB 
POSITION  VECTOR  FOR 
THIS  LEG 


i 


DETERMINE  \ 
:pa  FOR  THIS) 

leg  7 


no  (detection  does 
not  occur) 
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B 


yes  (detection  occurs) 


SHIP  IS  BEYOND  R 
FROM  PATROL  LINE  AT 
TIME  END  OF  THIS  LEG 
IS  REACHED 


0000  legs 
EXAMINED 


TAKE  ANOTHER 
SAMPLE 


<hK 


960  SAMPLES 
TAKEN 


7S 


no 


PRINT  RESULTS 
FOR  THIS  SPEED  RATIO 

CROSSING  POINT 
AVERAGE  NO.  LEGS 
SHIP  SPEED 
SPEED  RATIO 
RELATIVE  SPEEDS 
ROUTE  ANGLE 
PROB.  DETECTION  = 

(no.  detections) 
(no.  samples) 


1 


5 


SHIP  SPEED  =  90  KNOTS 


> 


no 


< 


i 


SHIP  SPEED  <  25  KNOTS 


o 


yes 


& 


INCREMENT  SHIP 
SPEED  BY  5 
KNOTS 


yes 
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BASIC  MODEL 


4 


INCREMENT 
SHIP  SPEED 
BY  10  KTS. 


,   V   N 
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EXTENDED  MODEL 


Q 


SHIP  SPEED  >  90  KNOTS 


* 


no 


no 


i 


4 


SHIP  SPEED  .   25  KNO 


yes 


4 


INCREMENT 
SHIP  SPEED 
BY  5  KNOTS 


yes 
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EXTENDED  MODEL 


Q 


SHIP  SPEED>  50  KNOT 


yes 


i 


5 


no 


CROSSING  POINT 
90  MILES  FROM 

CENTER  OF  SUB 
PATROL  LINE 


SET  SHIP 
SPEED  TO 
90  KNOTS 


4 


INCREMENT 
SHIP  SPEED 
BY  10  KNOTS 


no 


/  CR 


4 


CROSSING  POINT 


J  ^30  MILES  FROM 
CENTER  OF  SUB 
V    PATROL  LINE     J 


4 


INCREASE  CROSS 
POINT  DISTANCE 
FROM  CENTER  BY 
5  MILES 


£ 


INCREASE  CROSS 
POINT  DISTANCE 
FROM  CENTER  BY 
20  MILES 
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SUBROUTINE  WYEB 


no 


4 


DECREASE  SUB'S 
Y-COORDINATE 
BY  500  YARDS 


e 


& 


c 


ENTER 


3 


<J 


t 


SUB  HEADING 


"upy 


yes 


V 


INCREASE  SUB'S 
Y-COORDINATE 
BY  500  YARDS 


UB  REACHED  BOTTOM  \yes 
OF  PATROL  LINE 


no 


< 


yes/SUB  REACHED  TOP  OF 


PATROL  LINE 


REVERSE  SUB'S 
DIRECTION 


i> 


C 


no 


0 


EXIT 


") 
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SUBROUTINE  CPACHK 


C 


ENTER 


D 


COMPUTE  DOT  PRODUCT 
(REF  VECTOR) • (REL  SPEED  VECTOR) 


,     v     N 

f   DOT  PRODUCT  A 
V   NEGATIVE   ) 


no  (cpa  is  ref  vector 


yes  (cpa  not  reached  at 
start  of  leg) 


& 


COMPUTE  DOT  PRODUCT 
(TERMINAL  VECTOR) • (REL  SPEED  VECTOR) 


no 


I 


>•  \ 


DOT  PROD 
NEGATIVE 
OR  ZERO  J 
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i 


CPA  LIES  BETWEEN 
REFERENCE  VECTOR 
AND  TERMINAL  VECTOR: 
COMPUTE  CPA 


B 


V 


CPA  =  TERMINAL 
VECTOR 


^     EXIT     J 
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APPENDIX  C 

c  ************************* 

C  *  NUMERICAL  INTEGRATION  * 

C  *       BASIC  MODEL       * 

Q  ************************* 


CCMMCN/ONE/SGN,YB,ZN 
INTEGER  DALF 

DIMENSION  UV(2),W1V(2),W2V(2) , V1V( 2 > , V2V ( 2  )  f R1V( 2 ) , 
1RNV(2),UW1V(2)  ,UW2V(2),TEMV(2),V«V(2),VV(2),UWV(2) 

PI=3.1M5926536 

C  FORMAT  CONTRCL  CCUNTER  (INITIALIZE) 

IZ  =  2 

C  SUB  SPEED 

V  =  10. 

C  SAMPLE  SIZE 

JKL=960 

C  NP.  OF  ANGLE  (ALFA)  INCREMENTS 

IAI  =  9 

NA  =  5 

C  NO.  CF  RHC  INCREMENTS 

IBI=12 

C  DEFINE  BOTH  SUB  VECTORS 

CALL  VINS(V1V,0.,V) 
CALL  VINS(V2V,0.,-V) 

C  SIZE  OF  DETECTION  RADIUS  (CEFINITE  RANGE  LAW) 

R  =  3C. 

C  LENGTH  OF  SUB  PATRCL 

Pl=120. 

C  INITIAL  X-POSITICN  OF  PATRCL  LINE 

XP=R 

C  PRINT  FIRST  PAGE  HEADING 

WPITE(6,60C)R,PL,XP,V,JKL 

C  INITIAL  ANGLE 

ALF=0. 

C  LOOP  TO  INCREMENT  ALFA 

DC  8C0  IA=1,IAI 

C  CONVERSION  OF  RACIANS  TO  DEGREES 

IF(IA-1)105,23,21 
23  DALF=ALF 
GC  TC  22 

21  DALF=( IA+3)*5 

C  INITIAL  SHIP  SPEED 

22  U=10. 

C  CALCULATE  POSITION  OF  BCTTCM  CF  PATRCL  LINE  TO 

C  INSURE  CROSSING  AT  CENTER  OF  PATROL 

YP=XP*TAN( ALF)-(PL/2« ) 

C  LOOP  TO  INCREMENT  RHC 

DO  801  IB=1,IBI 
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C  INITIALIZE  LEG  CCUNT 

SLFG=0. 

CCSA=CCS(ALF> 
SINA  =  SIMALF) 

C  DEFINE  VECTOR  CF  SHIP  VELOCITY 

CALL  VINS(UV,U*COSA,U*SINA) 

C  COMPLTE  1/V 

VI=1./V 

C  COMPUTE  SHIP/SUB  SPEED  RATIO 

RHC=U*VI 

C  INITIALIZE  SUB  DIRECTIGN  FCR  DETERMINISTIC 

C  CALCULATION 

ZN=1. 

C  INITIALIZE  SUB  POSITION  FOR  DETERMINISTIC 

C  CALCULATION 

YB  =  YP 

C  INITIALIZE  DETECTION  COUNT 

DET=C. 

C  COMPUTE  MAXIMUM  TIME  FC*  CALCULATIONS 

TMAX  =  (XP«-P)/(U*COSA)-H. 

C  DETERMINE  RELATIVE  SPEED  VECTORS  AND  VALUES 

C  OF  MAGNITUDES  ANC  ANGLES 

CALL  VADD(VlVfUVfklV»-l) 

CALL  VADD(V2V,UVfW2V,-l) 

CALL  VVAL(W1V,W1,UW1V,BW1,111) 

CALL  VVAL<W2V,W2fUW2VtBW2,lll) 

C         DETERMINE  IF  DETECTION  OCCURS 

C  EACH  CCUNT  OF  THE  LOOP  INDEX  REPRESENTS 

C  ONE  SAMPLE 

DO  B  J=1,JKL 

C  DETERMINE  SUB  POSIT 

C  INPUTS 

C  XP:   X-DISTANCE  TO  PATROL  LINE 

C  YP:   Y-CCCRDINATE  OF  LOWER  LIMIT  CF  PATROL 

C  LINE 

C  PL:   LENGTH  OF  PATROL  LINE 

C  OUTPUTS 

C  PC:   RANGE  TO  SUB'S  INITIAL  POSIT 

C  BO:   BEARING  OF  SUB'S  INITIAL  POSIT 

C  YP  (VIA  CCMMCN/ONE/I :   INITIAL  Y-COORDINATE 

C  OF  THE  SUBMARINE 

C  R1V:  INITIAL  SUB  POSITION  VECTOR 
CALL  WYEB(XP,YP,PL,RC,BC,R1V) 

C  TIME  TO  ENC  CF  FIRST  (RANDCM)  LEG 

IF(SGN)7C0tl05,702 

C  TIME  TO  GO  FOR  SUB  TO  REACH  LIMIT  OF  PATROL  LINE 

C  IF  INITIAL  SUB  VELOCITY  IS  "DOWN" 

7CC  T1=(YB-YP)*VI 
GO  TO  703 

C  TIME  TO  GO  FOP  SUB  TO  REACH  UPPER  LIMIT  CF  PATROL 

C  LINE 

7C2  Tl=( YP+PL-YB)*VI 

7C3  T=Tl 
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C  INITIALIZE  TOTAL  TIME 

TN=0. 

C  INDICATE  THAT  THE  FIRST  (RANDOM)  LEG  IS  TO  BE  EX- 

C  AMINED 

MM=1 

C  COMPUTE  LENGTH  OF  RANDOM  LEG 

GC  TO  104 

C  LOOP  TO  INCREMENT  LEG  CCUNT 

103  DC  1C  IN=1,MM 

C  CHECK  TO  CETERMINE  INITIAL  DIRECTION  CF  SUB 

IFCSGNJ1C1, 105,102 

C  INITIAL  DIRECTION  IS  UP;  USE  PARAMETERS  FOR  THE 

C  UP  LEG 

101  Tfc=TK2 

CALL  C0PI(LWV,UW2V) 

W=W2 

CALL  C0PI(WV,V«2V) 

CALL  C0PI(VV,V2V) 

GC  TC  106 

C  INITIAL  DIRECTION  IS  DOWN;  USE  PARAMETERS  FOP  THE 

C  DOWN  LEG 

102  CALL    C0PKUWV,UK1V) 
TW=TW1 

W=W1 

CALL  C0PI(V*V,V<1V) 

CALL  C0PHVV,V1V) 

C  REVERSE  DIRECTION  CF  NEXT  LEG 

106  SGN=-SGN 

C  COMPUTE  RELATIVE  LENGTH  OF  THIS  LEG 

CALL  VMULT(UWV,TEMV,TM 

C  DETERMINE  VECTOR  FCR  FINAL  POSITION  ON  THIS  LEG 

CALL  VADD(P1V,TEMV,RNV,1) 
CALL  CPACHK(R1V,RNV,UV«V,RMIN,IN,  J) 

C  INCREMENT  LEG  CCLNT 

SLEG=SLEG*1, 

C  DETERMINE  IF  DETECTION  OCCURRED  ON  THIS  LEG 

IF(RMIN-P)3,3,11 

C  NO  DETECTION  THIS  LEG;  DETERMINE  IF  SHIP  IS  OUT 

C  OF  RANGE  BEYOND  PATROL  LINE 

11  TN=TN+T 
4  IF(TN-TMAX)6,8,8 

C  IF  THIS  IS  NCT  THE  FIRST  (RANDCMI  LEG,  CCNTINUE 

6  IF(MM-1)105,9,10 
10  CCNTINUE 
GO  TO  8 

C  LENGTH  CF  TIME  REQUIRED  FOR  SUB  TO  TRAVEL  ENTIRE 

C  SPAN  OF  PATROL  LINE 

9  T=PL*VI 

C  PROVIDE  FCR  POSSIBILITY  OF  10000  LEGS  MAXIMUM 

MM=10000 

C  COMPUTE  LENGTHS  OF  SECOND  AND  SUCCEECING  LEGS 

104  TW1=T*W1 
TW2=T*W2 
GO  TC  103 
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C  INCREMENT  CETECTICN  COUNTER 

3  CET  =  DET«-1. 
8  CONTINUE 
XJKL=JKL 

C  COMPUTE  PROBABILITY  OF  DETECTIGN 

P=DET/XJKL 
C  COMPUTE  AVERAGE  NO.  OF  LEGS  FOR  THIS  SCENARIO 

XLEG=SLFG/XJKL 

C  PRINT  RESULTS  FOR  EACH  SPEED  RATIO 

WRITE(6,6C1)XLEG,U,RHC,W1,W2,DALF,P 

C  INCREMENT  SHIP  SPEEC  BY  5  KNOTS  IF  RHC<7, 

C  THEN  INCREMENT  BY  10  KNOTS 

IF(IB-7)16,17,17 

16  U=U+5. 

GO  TO  801 

17  U=U*10. 
801  CCNTINUE 

C  FORMAT  OUTPUT  TO  INCLUDE  RESULTS  FOR  TWO 

C  VALUES  OF  ALFA  ON  FIRST  PAGE,  AND  THREE 

C  VALUES  OF  ALFA  ON  SUBSEQUENT  PAGES 

IF(IA-IZ)12, 13,105 

C  ONE  LINE  FEEC  BETWEEN  VALUES  OF  ALFA 

12  WRITE(6.602I 
GO  TO  20 

C  BEGIN  NEW  PAGE;  INSERT  COLUMN  HEADINGS 

13  IZ  =  IZ«-3 

IF( IA-9)19,800, 105 

19  WPITF<6,603) 

C  INCREMENT  ALFA  BY  25  DEGREES  ON  THE  FIRST 

C  ITERATION,  THEN  BY  5  DEGREES  ON  SUBSEQUENT 

C  ITERATIONS 

20  ALF=ALF*NA*(PI/36.) 

800  MA=1 

6CQ  FCRMAT(1H1.15(/),T16,'  PROBABILITY  OF  RANDOM  ENCOUNTER* 
1/T26, 'BASIC  M0DEL'///T6, 'DEFINITE  RANGE  LAW • /T8, • DETE • 
2,'CTION  RANGE:  ',F5.1,'  MILES '//T6, 'LENGTH  OF  SUB  PAT' 
3,'ROL  LINE:  ',F5.1,'  MILES' //T6 ,' X-DISTANCE  TO  PATROL' 
5,'  LINE:  SFS.l,'  MILES' //T6, 'SUBMARINE  SPEEC: ', F5. 1, 
5*  KN0TS'//T6,'NC.  OF  SAMPLES  TAKEN  FOR  EACH  SCENARIO:' 
6,'  • ,I4///1X,'AVG.N0.  SHIP  SPEED  RELATIVE  SPEEDS  • 
7, 'ROUTE  PROBABILITY  0F'/2X,'LEGS  SPEED  RATIO  UPLE' 
8,'G    OOWNLEG    ANGLE     DETECTION'//) 

601  F0RMAT(1X,F5.2,1X,2(2X,F4.1 )  ,2( 3X ,F5. 1 ) , 6X  ,  I  2 ,8X, F6.4 ) 

602  FORMAT(Z) 

6C3  FCRMAT(1H1,15(/).1X,'AVG.N0.  SHIP  SPEEC  RELATIVE  SP' 
1,'EEDS  ROUTE  PROBABILITY  0F'/2X, 'LEGS  SPEED  RATI' 
2,'G   UPLEG    DCWNLEG    ANGLE     DETECT ICN •// ) 

105  STOP 
END 
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************************* 

*  NUMERICAL  INTEGRATION  * 

*  EXTENDED   MODEL     * 
************************* 


CCMMON/ONE/SGN,YB,ZN 

DIMENSION  UV(2) » WLV< 2 ) «fe2V<2 ) .VI V < 2 ) , V2V ( 2 ) ,R1V < 2 ) , 
1RNV(2),UW1V<2> ,UW2V(2),TEMV(2),WV(2),VV<2),UWV(2) 

PI*3. 1415926536 

C  FORMAT  CONTRCL  COUNTER  (INITIALIZE} 

IZ=2 

C  SUB  SPEED 

V=10. 

C  SAMPLE  SIZE 

JKL=960 

C  NO.  CF  TRACK  CROSSING  INCREMENTS 

IAI=10 

C  NO.  OF  RHO  INCREMENTS 

IBI=11 

C  DEFINE  BOTH  SUB  VECTORS 

CALL  VINS(V1V,0.,V) 
CALL  VINS(V2V,0.»-VI 

C  SIZE  OF  DETECTION  RADIUS  (DEFINITE  RANGE  LAW) 

R*30. 

C  LENGTH  OF  SUB  PATROL 

PL=120. 

C  INITIAL  X-POSITICN  OF  PATROL  LINE 

XP=R 

C  PRINT  FIRST  PAGE  HEADING 

WRITE(6,6C0>R,PL,XP,V,JKL 

C  INITIAL  ANGLE 

ALF*0. 
YP*-PL/2. 

C  LOOP  TO  INCREMENT  CROSSING  POSITION 

DO  800  IA=1,IAI 

C  INITIAL  SHIP  SPEED 

U=10. 

C  CALCULATE  POSITION  OF  BOTTCM  OF  PATRCL  LINE  TO 

C  INSURE  CROSSING  AT  CENTER  OF  PATROL 

C  LOOP  TO  INCREMENT  RHO 

DO  801  IB=1,IBI 

C  DELETE  DATA  FOR  5<RH0<9 

IF(U-50. 720,20, 21 

C  GENERATE  OATA  FOR  RH0»9 

21  IF(U-90.)17, 22,105 

C  INITIALIZE  LEG  COUNT 

20  SLEG-O. 
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22  CCSA=CCS(ALF) 
SINA=SIN(ALF) 

C  DEFINE  VECTOR  CF  SHIP  VELOCITY 

CALL  VINS(UV,U*COSA,U*SINA) 

C  COMPUTE  1/V 

VI=1./V 

C  COMPUTE  SHIP/SUB  SPEEC  RATIO 

RHO=U*VI 

C  INITIALIZE  SLB  DIRECTION  FCR  DETERMINISTIC 

C  CALCULATION 

ZN  =  1. 
C  INITIALIZE  SUB  POSITION  FOR  DETERMINISTIC 

C  CALCULATION 

YP=YP 

C  INITIALIZE  DETECTICN  COUNT 

DET=0. 

C  COMPUTE  MAXIMUM  TIPE  FOR  CALCULATIONS 

TMAX  =  (XP*R)/(U*COSA)-H. 

C  CFTERMINE  RELATIVE  SPEED  VECTORS  AND  VALUES 

C  CF  MAGNITUDES  AND  ANGLES 

CALL  VADD(VlVfUV,WlVt-i) 

CALL  VADD(V2V,UV,V»2V,-1) 

CALL  VVAL(V.IV,  Wi ,UW1V  ,Bfcl  ,1 1 1  ) 

CALL  VVAL(W2V,W2,UW2V,BW2,11U 

C         DETERMINE  IF  DETECTICN  OCCURS 

C  EACH  COUNT  OF   INDEX  J  REPRESENTS  ONE  SAMPLE 

DO  8  J=1,JKL 

C  DETERMINE  SUB  POSITION 

C  INPUTS 

C  XP:   X-DISTANCE  TO  PATROL  LINE 

C  YP:   Y-CCCRDINATE  OF  LOWER  LIMIT  CF  PATROL 

C  LINE 

C  PL:   LENGTH  OF  PATROL  LINE 

C  OUTPUTS 

C  PC:   RANGE  TO  SUB'S  INITIAL  POSITION 

C  BO:   BEARING  CF  SUB'S  INITIAL  POSITION 

C  YB  (VIA  CCMMON/CNE/):   INITIAL  Y-COCRDINATE 

C  OF  SUBMARINE 

C  P1V:  INITIAL  SUB  POSITION  VECTOR 
CALL  ViYEB(XP,YPtPL,RC,BC,RlV) 

C  TIME  TO  END  OF  FIRST  (RANDCMI  LEG 

IF(SGN)70C, 105,702 

C  TIME  TO  GO  FOR  SUB  TO  REACH  LIMIT  OF  PATROL  LINE 

C  IF  INITIAL  SUB  VELOCITY  IS  "DOWN" 

7CC  Tl=( YB-YP)*VI 
GO  TO  703 

C  TIME  TO  GO  FOP  SUB  TO  REACH  UPPER  LIMIT  CF  PATROL 

C  LINE 

7C2  T1=(YP+PL-YB)*VI 

703  T=T1 

C  INITIALIZE  TOTAL  TIME 

TN=0. 

C  INDICATE  THAT  THE  FIRST  (RANDOM)  LEG  IS  TO  BE  EX- 

C  AVINED 

MM*1 
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C  COMPUTE  LENGTH  OF  RANDOM  LEG 

GO  TO  104 

C  LOOP  TO  INCREMENT  LEG  COUNT 

103  DO  10  IN=1,MM 

C  CHECK  TO  DETERMINE  INITIAL  DIRECTION  CF  SUB 

IF<SGN)101, 105,102 

C  INITIAL  OIRECTION  IS  UP;  USE  PARAMETERS  FOR  THE 

C  UP  LEG 

101  TW=TW2 

CALL  C0PI(UWV,UV«2V) 

W-W2 

CALL  C0PI(WV,W2V) 

CALL  C0PI(VV,V2V) 

GO  TO  106 

C  INITIAL  DIRECTION  IS  DOWN;  USE  PARAMETERS  FOR  THE 

C  DOWN  LEG 

102  CALL    COPHUWVvUMlV) 
TW^TWl 

W=W1 

CALL  COPI(WV,klV) 

CALL  C0PI(VV,V1V) 

C  REVERSE  DIRECTION  CF  NEXT  LEG 

106  SGN=-SGN 

C  COMPUTE  RELATIVE  LENGTH  OF  THIS  LEG 

CALL  VMULT(UWV,TEMV,TW» 

C  DETERMINE  VECTOR  FOR  FINAL  POSITION  ON  THIS  LEG 

CALL  VADD(R1V,TEMV,RNV,1) 
CALL  CPACHK(RlV,RNV,UMV,RMINvINvJ) 

C  INCREMENT  LEG  CCUNT 

SLEG«SLEG+1. 

C  DETERMINE  IF  DETECTION  OCCURRED  ON  THIS  LEG 

IF<RMIN-RJ3,3,11 

C  NO  DETECTION  THIS  LEG;  DETERMINE  IF  SHIP  IS  OUT 

C  OF  RANGE  BEYOND  PATROL  LINE 

11  TN*TN+T 
4  IF<TN-TMAX)6,8,8 

C  IF  THIS  IS  NOT  THE  FIRST  (RANDOMI  LEG,  CONTINUE 

6  IF<MM-1)105,9,1G 
10  CONTINUE 
GO  TO  8 

C  LENGTH  OF  TIME  REQUIRED  FOR  SUB  TO  TRAVEL  ENTIRE 

C  SPAN  OF  PATROL  LINE 

9  T»PL*VI 

C  PROVIDE  FOR  POSSIBILITY  OF  10000  LEGS  MAXIMUM 

MH-10000 

C  COMPUTE  LENGTHS  OF  SECOND  AND  SUCCEECING  LEGS 

104  TW1=T*W1 
TW2»T*W2 
GO  TO  103 

C  INCREMENT  DETECTION  COUNTER 

3  DET-DET+1. 
8  CONTINUE 
XJKL«JKL 
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C  COMPUTE  PRCBABILITY  CF  OETECTIQN 

P=DET/XJKL 

C  COMPUTE  AVERAGE  NO.  CF  LEGS  FOR  THIS  SCENARIO 

XLEG=SLEG/XJKL 

C  GENERATE  CROSSING  COORDINATE  FCR  PRINTOUT 

XYP=-YP 

C  PRINT  RESULTS  FOR  EACH  SPEED  RATIO 

WRITE (6, 60 II XLFG,U,RHC,W1,W2,XYP,P 

C  INCREMENT  SHIP  SPEEC  BY  5  KNCTS  IF  RHCO, 

C  THFN  BY  10  KNCTS 

IF(IB-5)16,17,17 

16  U  =  U«-5. 

GC  TC  801 

17  U=U+1C. 
8C1  CONTINUE 

C  FORMAT  OUTPUT  TO  INCLUDE  RESULTS  FCR  TWC 

C  VALUES  OF  ALFA  ON  EACH  PAGE 

IF(IA-IZ)12tl3,105 

C  ONE  LINE  FEEC  BETWEEN  VALUES  OF  ALFA 

12  WPITE<6,602) 
GO  TC  18 

C  BEGIN  NEW  PAGE;  INSERT  COLUMN  HEADINGS 

13  IZ=IZ+2 
IF(IA-10)19. 18,105 

19  WPITE(6,6C3) 


C 
C 
C 


INCREMENT  CROSSING  POSITION  BY  5  MILES 
UNTIL  WITHIN  30  MILES  OF  PATROL  TERMINUS 
THEN  INCREMENT  BY  20  MILES 
19  IF(IA-7)14,15, 15 

14  YP=YP-5. 
GO   TO    800 

15  YP=YP-20. 

\n.r>    rnMTTMiic 


800  CONTINUE 


6CC  FCRMATdHl,  8(  /  )  ,T16  ♦  •  PRCBABI  LI  TY  CF  RANCCM  ENCOUNTER* 
1/T24, 'EXTENDED  KDOEL' 

1///T6, 'DEFINITE  RANGE  LAW '/T8, • DETECT  ION  RANGE:  '.F5.1 
2,'  MILESV/T6, 'LENGTH  CF  SUB  PATROL  LINE:  ',F5.1,»  MI" 
3,  'LES'//T6,'X-DISTANCE  TO  PATROL  LINE:  ,tF5.lf«  MILES' 
4//T6, 'CROSSING  ANGLE:  CO  OEGREES  (PERPENDICULAR)' 
5//T6, 'SUBMARINE  SPEED  :•  ,F5. 1 ,  •  KNCTS • //T6, • NO.  OF  SAM' 
6,'PLES  TAKEN  FOR  EACH  SCENARIO:  • , I4///1X , • AVG.NO.  SH' 
7, 'IP  SPEED  RELATIVE  SPEEDS  CROSS  PROBABILITY  OF'/ 
82X,'LEGS  SPEED  RATIC  UPLEG  DOWNLEG  POSIT  DE' 
9,'TECTION'//) 

601  FCRMAT(1X,F5.2,1X,2(2X,F4.1  ) , 2( 3X, F 5. 1 ) , 5X, F5. I, 7X, 
1F6.4) 

602  FORMAT!/) 

503  F0PMAT<1H1.17(/I.IX1'AVG.N0.  SHIP  SPEED  RELATIVE  SP» 
1,'EEDS  CROSS  PROBABILITY  0F'/2X,'LEGS  SPEED  RATI' 
2,'C   UPLEG    DCWNLEG    POSIT     DETECTION'//) 

105  STOP 
END 
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c  ************** ****** ************ 

C  *  NUMERICAL  INTEGRATION  MOCELS  * 

C  *  SUBROUTINES  * 

£  ******************************** 


SUBPOLTINE  WYEB<XP,YP,PL,RC,8C,A) 

C  THIS  SUBROUTINE  GENERATES  SECUENTIAL  V- 

C  COORDINATES  ON  A  SPECIFIED  LINE  SEGMENT; 

C  fcHEN  AN  ENC  PCINT  OF  THE  LINE  IS  GENERATED, 

C  THE  DIRECTION  OF  THE  SEQUENCE  IS  REVERSEC 

C  INPUTS 

C  YB=Y-POSITIGN  CF  SUB  ON  PATROL  LINE 

C  SGN=INITIAL  CIRECTION  OF  SUB 

CCMMON/ONE/SGN,YB,ZN 
DIMENSION  A(2) ,E(2) 

C  INCREMENT  SUB  POSITION  BY  500  YDS  IF  ZN=1; 

C  DECREMENT  SUB  POSIT  BY  500  YDS  IF  ZN  =  -1 

YB=YB+0,25*ZN 

C  REVERSE  SUB  COURSE  WHEN  IT  REACHES  ENC  OF 

C  PATROL  LINE 

IF(ZN*(2,*(YP-Y8)+PL)4PL)9,9,10 

C  POSITION  SUB  EXACTLY  AT  TERMINUS  OF  PATROL 

C  LINF  PRIOR  TC  TURN 

9  YB=YP+0.5*FL*(1.+ZN) 

C  REVERSE  CCURSE  CF  SUB 

ZN=-ZK 
10  SGN=ZN 

C  CETEPMINE  INITIAL  RANGE  AND  BEARING  TC  SUB 

A(1)=XP 
A<2)=YB 

CALL  VVAL(A,R0,B,8C,101) 
RFTUPN 
END 
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SLBRCLiTINE    VINS(A,X,Y) 

C  THIS  SUBROUTINE  INSERTS  THE  X  ANC  Y  VECTOR 

C  COMPONENTS  INTC  ARRAY  A 

DIMENSION  A<1) 

A(1)=X 

A(2)=Y 

RETURN 
END 


SUBROUTINE  COPKA.BI 

C  THIS  SUBROUTINE  GENERATES  THE  VECTOR  A  SUCH  THAT 

C  FOR  A  GIVEN  VECTOR  B,  A  =  B 

DIMENSION  A(l) ,P(1) 
DC  1  1=1.2 
1  A(I)=B(I) 

RFTUPN 
END 


SUBROUTINE  VADD( A ,B ,R ,N) 

C  THIS  SUBROUTINE  PERFORMS  THE  FOLLOWING  VECTCR 

C  OPERATIONS: 

C  IF  N>C:  R=A*>6  (VECTOR  ADDITION) 

C  IF  N<C:  R*A-8  (VECTOR  SUBTRACTION) 

C  IF  N=0:  THE  DOT  PRODUCT  OF  A  ANO  B  IS  FCRMEC  AND 

C  RESULT  IS  PLACED  IN  THE  FIRST  MEMBER  CF  ARRAY  R 

DIMENSION  A(l),e(l),R(l) 
IF(N)l,2,3 

C  ADD  VECTORS 

3  DO  4  1*1.2 

4  R  ( I )=A( I »*B( I ) 

RETURN 

C  SUBTRACT  VECTCRS 

1  DC  5  1=1.2 

5  R(I)*A(I)-B(I) 

RETURN 

C  COMPUTE  DCT  PRODUCT  OF  VECTORS 

2  DO  6  1=1,2 

6  R(I)=AU)*B(I) 
R(1)=P<1)*P<2) 

RFTUPN 
END 
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8 


SUBROUTINE  VMULT<A,B,X> 

THIS  SUBROUTINE  MULTIPLIES  THE  VECTOR  A  BY  THE 
SCALAR  X  TC  FORM  THE  VECTOR  B 

DIMENSION  A(l),e(l) 
DO  1  1=1.2 

BU)=X*A(I) 

RETURN 

END 


SUBROUTINE  VVAL( A f VL,UV, ANGLE  .M ) 

C  THIS  SUBROUTINE  DETERMINES  THE  FOLLOWING 

C  VALUES  OF  A: 

C  A=INPUT  VECTOR  (SUPPLIED  BY  THE  CALLING  PROGRAM) 

C  VL*MAGNITUDE  OF  A 

C  UV=UNIT  VECTOR  PARALLEL  TO  A 

C  ANGLE=ANGLE  OF  A  MEASURED  CCW  FROM  X-AXIS 

DIMENSION  Am»UV(l> 
IF(M-100)2,3,4 

4  !F<M-110)5t6,7 

C  COMPUTE  VECTOR  LENGTH  ONLY 

3  VL=SQRT<A(l»**2*At2)**2l 
RETURN 

C  COMPUTE  ANGLE  OF  VECTOR  ONLY 

2  ANGLE=ATAN2<A<21,M1)  ) 
RETURN 

C  COMPUTE  BOTH  LENGTH  AND  ANGLE  OF  VECTOR 

5  VL»SQRT(A(1!**24A(2)**2> 
GO  TO  2 

C  GENERATE  UNIT  VECTOR  ONLY 

6  VL«SQRT<A<1I*«2*A<2)**2) 
DO  1  1=1,2 

1  UVm=A(I)/VL 
RETURN 

C  COMPUTE  ANGLE  AND  GENERATE  UNIT  VECTOR 

7  ANGLE=ATAN2(A(2)(A(1H 
GO  TO  6 

END 
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SUBROUTINE  CPACFK  {  R  1 ,  PN  , UWV  ,  RCPA  ,  IN,  J  ) 

C  THIS  SUBROUTINE  DETERMINES  THE  CLOSEST  POINT 

C  TO  THE  ORIGIN  OF  ANY  LINE  SEGMENT.  THE  END 

C  POINTS  OF  WHICH  ARE  DEFINEC  BY  VECTORS 

C  R1=RFFEPENCE  VECTOR 

C  RN=TFRMINAL  VECTOR 

C  UWV=UMT    VECTOR    PARALLEL    TO    CRM 

C  RCPA=DISTANCE  AT  CLOSEST  POINT  OF  APPROACH  (CPA) 

C  IN=LEG  NO.  BEING  EXAMINED  AT  THIS  EXECUTION 

C  J=SAMPLF  NO.  BEING  EXAMINED  AT  THIS  EXECUTION 

DIMENSION  FL(1) tRN(l) ,UWV<1) ,X(2),Y(2),CPA(2),Z<2), 
1PT(2) 
CCMMCN/TWC/Y,N 

C  COMPLTE  DCT  PRCCLCT:  R1*UWV=X 

CALL  VADD(PI,UWV,X,0) 
IF(X(l))i,2,2 

C  CPA  WILL  CCCUR  ON  THIS  LEG;  COMPUTE  THE  COT  PROD- 

C  LCT:  RN*UWV=X 

1  CALL  VADD(RN,UWV,X,0) 
IF(X<1))3,3,4 

C  CPA  WILL  NOT  CCCUR  ON  THIS  LEG;  CPA=R1 

2  CALL  VVAL(P1, RCPA, X,C, 100) 
RFTURN 

C  CPA  WAS  NOT  REACHED  ON  THIS  LEG;   CPA=RN 

3  CALL  WAL(PN, RCPA, X,D, 100) 

RFTURN 

C  CETEPMINE  CPA;  ASSUME  LEG  HAS  LENGTH  <  250  MILES 

4  DO  5  1=1, 5C0 

C  USE  2-MILE  INCREMENTS 

P=2.*I 

C  INCREMENT  CRM  BY  2  MILES 

CALL  VMULT(UWV,Y,P) 

C  ADD  TO  REFERENCE  VECTOR 

CALL  VADD(P1,Y,X,1) 

C  COMPUTE  DCT  PRCCUCT   X*UWV=Z 

CALL  VADD(X,UWV,Z,0) 

C  IF  Z  BECOMES  POSITIVE,  MAKE  FINE  INCREMENT  CHECK 

C  BELOW;  IF  Z(1)=0,  X=CPA;  OTHERWISE  TRY  AGAIN 

IF ( Z ( 1 ) )5,6,7 

C  INCREMENT  CRM  BY  0.1  MILE 

C  PRESERVE  PRECECING  VECTOR 

7  CALL  VMULT(UWV,Y,P-2.  » 
CALL  VADC(P1,Y,X,1) 

C  CHECK  LAST  2  MILES  OF  DRM  IN  O.l-MILE  INCREMENTS 

DC  8  1=1,21 
P=0.1*I 

C  INCREMENT  CRM  BY  0.1  MILE  AND  FIND  CPA 

CALL  VMULT(UWV,Y,P) 

C*LL  VADD(X,Y,RT,1) 

CALL  VADD(RT ,UWV,Z,0) 
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IF(Ztl))8,6,6 

8  CONTINUE 

C  ERROR  EXISTS  IN  PROGRAM  IF  NEXT  INSTRUCTION  IS 

C  EXECUTED 

WRITE(6,61)IN,J 
61  F0RMAT(//10X,»LEG«.1X,I4,1X, • CF  SAMPLE* , 15, IX, • GENER* , 
l'ATED  ERROR  IN  CPACHK  ROUTINE1//) 
IF(Z(1))  5,6,7 

5  CONTINUE 

C  ORM  LEG>250  MILES  IF  NEXT  INSTRUCTION  IS 

C  EXECUTED 

WRITE(6,60)INtJ 
60  F0RMAT(//10X,*LEG«.1X1I4,1X,»0F  SAMPLE* • 15, IX, • IS 
1GREATER  THAN  250  MILES'//) 
RETURN 

C  COMPUTE  CPA  DISTANCE 

6  CALL  VVAL(X,RCPA,Y,D,100) 
RETURN 

END 
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APPENDIX  D 

f  *********************** 

f  *    CAMPCM  SAMPLING    * 

*   MIDPOINT  CROSSING   * 
p  ***  ************** ****  ** 

Tf  MMPNi/n^F/SGNfYB 

DIMENSION  liV(2)fWlV<2)»W2V(2»tVlV(2» , V2V( 2  )  t R  1  VI  2  )  , 

ieNvm,uwiv<2)  fuw2v<2 )  f  TEMvm.wvm  tvvm  tuwv(2) 

C  iriTI/l  T^F    URN    GENFRATOR 

pr*STT=UPN(C) 

PT=3,  1  -lc92^e;3^ 
C  MIR    SPFFn 

C  SAMPLE    CI7F 

JK! -icon 

f  Mr.     OF    *NGLF     (ALFA)     INCREMENTS 

IAI=17 

C  NO.     hf     FHH    INCRFMFNTS 

TPI  =  o 

C         DEFINE  *OTH  SUP  VcCTOPS 

C'LL  VINS(VlVf0.  »V) 
CALL  VTNS(V2V,0.,-V) 

f         SIZr  OF  HETPCTION  PAPIUS  (DEFINITE  RANGE  LAW) 

P=3C. 
C         LENGTH  TF  SUP  PATPHL 

Pl=120. 
C         INITIAt  X-PnsiTION  OF  patrol  lin: 

Xf=P 

WFITFC^tfOr)R,PL,XPtV,JKL 
r  IMTIAl  ANGl  c 

A.I  p  =  C.^ 
f.         t  TP  TO  INCREMENT  ALFA 

or  poo  U  =  l  f  tai 
r         INITIAL  SHIP  SPEED 

C         CALCULATE  PHSITION  OF  POTTO*  of  PATR°L  LINE  to 
r  INSURF  CROSSING  AT  CENTER  OF  PATROL 

YF=Xp*TAN< ALF)-(PL/2. ) 
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C         LOOP  TO  INCREMENT  RHO 

DO  801  IB=1,IBI 
C         INITIALIZE  LEG  COUNT 

SLEG=0. 

COSA=COS(ALF) 
SINA=SIN(ALF) 

C         DEFINE  VECTOR  OF  SHIP  VELOCITY 

CALL  VINS(UV,U*COSA,U*SINA) 
C         COMPUTE  SHIP/SUB  SPEED  RATIO 

RHO=U/V  { 

C         INITIALIZE  SUB  DIRECTION  FOR  DETERMINISTIC 
C         CALCULATION 

SGN-1. 
C         INITIALIZE  DETECTION  COUNT 

DET-O. 
C         COMPUTE  MAXIMUM  TIME  FOR  CALCULATIONS 

TMAX  =  (XP+R)/(U*COSA)  +  U 

C         DETERMINE  RELATIVE  SPEED  VECTORS  AND  VALUES  OF 
C         MAGNITUDES  AND  ANGLES 

CALL  VADD(VlV,UV,WlVt-l) 
CALL  VADD(V2V,UV,W2V,-l) 
CALL  VVAL(WlV,Wl,UVaV,BWl,lll) 
CALL  VVAL(W2V,W2fUW2V,BW2,lll) 

C         DETERMINE  IF  DETECTION  OCCURS 

C         EACH  COUNT  OF  THE  LOOP  INDEX  REPRESENTS  ONE  SAMPLE 

DO  8  J=1,JKL 

C  DETERMINE  SUB  POSIT  (UNIFORMLY  DISTRIBUTED) 

C  INPUTS 

C  XP:  X-DISTANCE  TO  PATROL  LINE 

C  YP:  Y-COORDINATE  OP  LOWER  LIMIT  OF  PATROL 

C  LINE 

C  PL:  LENGTH  OF  PATROL  LINE 

C  OUTPUTS 

C  RO:   RANGE  TO  SUB'S  INITIAL  POSIT 

C  BO:   BEARING  OF  SUB'S  INITIAL  POSIT 

C  YB  (VIA  COMMON/ONE/I :   INITIAL  Y-COORDINATE 

C  SUB 

C  R1V:  INITIAL  SUB  POSITION  VECTOR 

CALL  WYEB(XPfYP,PL,RO,BO,RlV) 

C         TIME  TO  END  OF  FIRST  (RANDOM)  LEG 

IF(SGN)700, 105,702 

C         TIME  TO  GO  FOR  SUB  TO  REACH  LIMIT  OF  PATROL  LINE 
C         IF  INITIAL  SUB  VELOCITY  IS  "DOWN" 

700  T1=(YB-YP)/V 
GO  TO  703 

C         TIME  TO  GO  FOR  SUB  TO  REACH  UPPER  LIMIT  OF  PATROL 
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r  |  Tfr 

?r?  T] =( yp+FI-YR) /V 

■*  c  i  T = T  i 
c       initialize  total  ttme 

tn*c« 
f       imoic^tt  that  the  first  (random)  leg  is  to  "f 

r  FXfrMINET 

C         COMPUTE  LENGTH  OF  RANDOM  LEG 

C         LOOP  TP  INCREMENT  LEG  COUNT 

Iq-2   rr   ic    in=i,mm 

C         CVfC*    TC  DETERMINF  INITML  DIRECTION  OF  SUR 

TF( SGNI1C1 ,105,102 

C         IN'ITjAi  OIRFCTION  IS  UP;  USC  PARAMETERS  «=0R  THE 
C         UF  LFG 

101  TW=TW? 

CfLL    COPI (UWV,UW2V) 
W  =  W2 

CALl  Cnpi(wv,W2V) 
CALL  C0PI(VV,V2V) 
GT  TP  106 

r         INITIAL  DIRECTION  TS  DOWN;  USE  PARAMETERS  FDR  THE 
C         DOWN  IFG 

102  CALl  CnPI (UUVfUWlVI 
TV-  =  TW1 

W  =  N1 

CALl  CnpI(WV,WlV) 

CALL  C0PI(W,V1V) 

C         PFVPPSE  DIRECTION  CF  NEXT  LCG 

10*  SFN  =  -SGN 

C         COMPUTF  RELATIVE  LENGTH  0«=  THIS  LEG 

CtLl     VMULT(UWV,TF*V,TW| 

C         DETERMINE  VECTOR  ECR  FINAL  PPSITION  ON  THIS  LEG 

CALl  VAD0(P1V,TEMV,RNV, H 

CALL  CPACHK(PlV,RNV,UWV,RMIN,IN, J) 

INCREMENT  LEG  COUNT 

SLEG*SLEG*1. 

C  DETERMINE    IP    DETECTION    OCCURRED    ON    THIS    LEG 

IF(PMIN-P)3f3fH 

11  tn*tn+t 

a    TF(TN-TMAXI6,<?,8 
*    IF(MM-j  nC5,4f  10 

1  1    CFNTINUF 
DC     Tf    R 
->    T  =  RL/V 
MM*10000 
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104  TW1=T*W1 
TW2=T*W2 
GO  TO  103 

3  DET=0ET+1. 
8  CONTINUE 

XJKL=JKL 

P=DET/XJKl 

XLEG=SLFG/XJKL 

WPITE(6f601)XLEG,U,RHC,Wl,W2,ALF,P 
801  U=U+10. 

WPITE(6f602) 
602  FORMAT(Z) 
800  ALF=ALF*PI/36. 

600  F0RMATUH1/////T161'PR0BABILITY  OF  RANDOM  ENCOUNTER1,/ 
1//T6, 'DEFINITE  RANGE  LAW /T8, 'DETECTION  RANGE:  NES.!. 
2'  MILES'//T6,  'LENGTH  OF  SUB  PATROL  LINE:  '.F5.1.'  MIL' 
3t'ES'//T6,'X-0ISTANCE  TO  PATROL  LINE: • tF5.it •  MlLES'// 
4T6,'SUBMAPINE  SPEED: '.F5. 1. •  KNOTS' //T6, 'NO.  OF  SAMPL • 
5,'ES  TAKEN  FOR  EACH  SCENARIO:  'fI4///T3,'  AVG.  NO.  SH' 
6, 'IP   SPEED  • 

7, 'RELATIVE  SPEED   ALFA    PROB.  DETECTION' /T5, 'LEGS    • 
8t«  SPEED  RATIO  UPLEG   DOWNLEG'//) 

601  F0RMAT(T4fF5.2t2X,3(lX,F5.1),2X,F5.1,3X,F6.4,7X,F6.4> 

105  STOP 
END 


SUBROUTINE  WYEB(XP,YP,PL,ROt BO, A) 

C         INPUTS 

C  YB=Y-POSITION  OF  SUB  ON  PATROL  LINE 

YB  IS  UNIFORMLY  DISTRIBUTED 

SGN=INITIAL  DIRECTION  (EQUALLY  LIKELY)  OF  SUB 

VELOCITY 

C0MM0N/0NE/SGN,YB 
DIMENSION  A(2)9B(2I 

DETERMINE  SUB  POSIT  RANDOMLY 

YB*YP+URN(1)*PL 

DETFRMINE  INITIAL  DIRECTION 


4 
6 

IF(URNU)-0.5)5t6,7 
GO  TO  4 

UP 

5 

SGN=1. 
GO  TO  8 

DOWN 

7 

SGN*-1. 

DETERMINE  INITIAL 

RANGE  AND  BEARING  TO  SUB 


8  A(1)=XP 
A(2)«YB 
CALL  VVAL(A,R0,B,B0,101> 

RETURN 
END 
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SURPPUTINf  CPACHK(RlfRN,UWV,RCPAf IN, J) 

riMEKSION  RlllliRNIll ,UWV(1 ),X(2),Y(2),CPAm,?(2) 
rrMMpK/TWO/V,N 

c       F  i=repfpence  VECTOR 

C         P>  =  TcpminAL  VECTOR 

r         UkV=UNIT  VECTOR  CF  ppm 

CClMPUTi     HOT    PRODUCT:    R1*IJWV=X 

rai    vapp(ri,uwv,x,o) 

IMX(  1)  )  UP, 2 

f         CPA  WTU  OCCUR  ON  THIS  LEO:  COMPUTE  THE  TOT 
f         PRODUCT:  RN*UWV=X 

1  C^Ll  VADP(PN,UWV,X,0) 
IMX(l)  13,3.4 

C         CPA  will  NOT  OCCUR  ON  THIS  LFG:  CPA=P1   cr»w°UTF 
r         CFA 

?  call  vv/l(pi ,rcpa,x,p,ioo) 

PFTUFK 
C         CFA  WAS  NHT  RFACHFT  DN  THI*  LEO:   CPA  =  RN:  COMPUTE 

r       c  p  a 

3  '"AH  VVAl  (FN,RCPA,X,D,100) 
BETUPN 

r         DETERMINE  CPA:  ASSUME  LEO  HAS  LENGTH<  250  MILES 

4  op  F  1=1, Fro 

C         UCF  HALF-HUF  INCRFMFNTS 

P=I/2. 
C         INCREMENT  DPM  R-Y  1/2  "ILE 

C«H  VMUtT(UWV,Y,P) 
C         £PP  TO  PEFPRENCE  VFCTOR 

CALL  VAPP(P1 ,Y,X,l) 
C         COMPUTF  DOT  PRODUCT   X  DOT  UWV^Z 

CALL  VAPP(X,UWV,Z,0) 

r         IF  7  CHANGES  SIGN.  OR  RECQMFS  ZFRO  VECTOR  X=CPA: 
F         PTHFPWIS'!  TRY  AGAIN 

T  F  (  7  C  I )  I  *,<>,* 

ff  CONTINUE 

C  PPM   LFG>250    MILFS 

WPITF(6,6C) IN, J 
6*    FPPMAT(//KX,  'LEG'  .  1 X.I4, IX , • OF    SAMPLE  •  ,  I  5  , 1  X  ,  •  I  S 
1GPE/TFF    TH«N    2«50    m|lF$«//) 
RFTUPM 
a    C*LL    VVAl  (X,RCPA,Y, 0,100) 
PFTUPM 
END 
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SUBROUTINE  VINS(AfX,Y) 

C         THIS  SUBROUTINE  INSERTS  THE  X  AND  Y  VECTOR 
C         COMPONENTS  INTO  ARRAY  A 

DIMENSION  A<1) 

A(1)=X 

A(2)=Y 

RETURN 
END 


SUBROUTINE  COPKA.BI 

C         THIS  SUBROUTINE  GENERATES  THE  VECTOR  A  SUCH  THAT 
C         FOR  A  GIVEN  VECTOR  B,  A*B 

DIMENSION  A(l), 8(1) 
DO  1  1*1.2 
1  A(I)*BUI 
RETURN 

END 


SUBROUTINE  VADD( A , B, R , N) 

C  THIS  SUBROUTINE  DOES  THE  FOLLOWING  VECTOR 

C  OPERATIONS: 

C  IF  N>0:  R=A+B 

C  IF  N<0:  R=A-B 

C  IF  N=0:  THE  OOT  PRODUCT  OF  A  AND  B  IS  FORMED  AND 

C  RESULT  IS  PLACED  IN  THE  FIRST  MEMBER  OF  ARRAY  R 

DIMENSION  All) «B(1)*RC1I 
IF(N)1,2,3 

C         ADD  VECTORS 

3  DO  4  1=1,2 

4  RU)«A(I!*Bm 

RETURN 

C         SUBTRACT  VECTORS 

1  DO  5  1*1,2 

5  R(I)=A(I)-B(I) 
RETURN 

C         COMPUTE  DOT  PRODUCT  OF  VECTORS 

2  DO  6  1=1,2 

6  R(I)=A(I)*B(I) 
R(1)=R(1)*R(2) 

RETURN 

ENO 


105 


SIJPPriJTINF     VMULT(4,R,X) 

THIC     ?«|PPOi|TINF    MULTIPLIES    THF    VECTOR     A    RY    the 
SCtUP     y    TO    FORM    THF    VECTOR    R 

DIMENSION    A(  1  )  ,B(  1) 
DP    1     1=1.2 

i  pm  =  x*Mn 

crTUPN 


EUFROUTINF  VVAL(  Af  VLf  UVtANGLF.  fM) 

THIS  5UPRPUTINE  DETERMINES  THE  FOLLOWING 
V£L»»FC  PP  A: 

A=INPUT  VECTOR  (SUPPLIED  FY  TH^  CALLING  PROGRAM) 

VL=M*GNITUOF    OF    A 

UV=UMT    VECTOR    PARALLEL    TO    A 

*NGLC=ANGIE    OF    A    MEASURED    CCW    FROM    X-AX*S 

D I  MENS  I  CM  All)  fUVU! 

TF(w-i00)2,3,4 
u    IF<M-llC)?t6f7 

COMPUTE  VECTOR  LENGTH  ONLY 

"*  Vl-SOPT(  Am**2+A(2)**2) 

PFTUPM 

COMPUTE  ANGLE  OF  VECTOR  ONLY 

2    ANGLF  =  ATAN2( A(2) ,  A(LM 
PFTIIPN 

COMPUTE  ROTH  LENGTH  AND  ANGLE  OF  VECTOR 

r    VL=SOPT( A( 1)**2*A(2)**2) 
or   to  2 

GENERATF  UNIT  VECTOR  ONLY 

6  VL=SOPT(A(l)**2+A(2)**2) 

pr   i  l*l .2 

1  UV(T )=A(t)/VL 

FETUFN 

COMPUTF  ANGLF  AND  GENERATE  UNIT  VECTOR 

i    £MGIE»ATA*2< A(2),A(11  ) 

GO  TP  h 

EMr 
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c  ***************************** 

C  *       RANDOM  SAMPLING       * 

C  *   NORMALIZED  CONVOY  ROUTE   * 

C  ***************************** 


COMMON/ONE/SGNtYB 

1V(2).W2V(2).V1V(2).V2. 

(2),UWV(2) 


DIMENSION  UV(2),W1V(2),W2V(2),V1V(2),V2V(2),RIV(2>, 
lRNV(2»,UWlV(2),UW2V(2l,TEMV(2),WV(2),VVr- 


REAL  NORM 

INITIALIZE  URN  GENERATOR 
POSIT=UPN(0) 
PI=3. 1415926536 

SUB  SPEED 
V=10. 

SAMPLE  SIZE 
JKL-10OO 

NO.  OF  ANGLE  (ALFA)  INCREMENTS 
IAI=17 

NO.  OF  RHO  INCREMENTS 
IBI=9 

DEFINE  BOTH  SUB  VECTORS 

CALL  VINS(V1V,0.,V) 
CALL  VINS(V2V,0.,-V) 

SIZE  OF  DETECTION  RADIUS  (DEFINITE  RANGE  LAW) 
R=30. 

LENGTH  OF  SUB  PATROL 
PL=120. 

INITIAL  X-POSITION  OF  PATROL  LINE 
XP*R 
WRITE(6,600)R,PLfXP,VfJKL 

INITIAL  ANGLE 
ALF=0.0 

LOOP  TO  INCREMENT  ALFA 
DO  800  IA*ltIAI 

INITIAL  SHIP  SPEEO 
U=10. 

LOOP  TO  INCREMENT  RHO 
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DP  PCI     !P  =  l,m 
C.         INITIALIZE  LEG  COUNT 

SL^c-D. 

rrcA=rn<(  AL^) 
SINA=SIN( ALF) 

C         DEFINE  VECTOR  OF  SHIP  VELOCITY 

CfLl  VT>,c(UV,U*CnSA,U*SINA) 

r  COMntiTF    SHIP/SUR    SPEED    PATm 

P»->p=U/V 

C  INITIALIZE  SUB  DIRECTION  FOP  DETERMINISTIC 

C       mcui£T!nN 

SGN=1. 
r  INITIALIZE  DETECTION  COUNT 

nfcT.c. 

C         COMPUTE  MAXIMUM  TIME  FOR  CALCULATIONS 

TVAX  =  (XP  +  P|  /(U*COSA)+l. 

C         DFTFPVINE  RELATIVE  SPEED  VECTORS  AND  VALUES  HF 
C         MAGNITUDES  AND  ANGLES 

CALL  VADP(V1V,UV,W1V,-1) 

CALl  VADD(V2V,UVfW2V,-l) 

C^LL  VVAL(W1V, WltUWlVtBWl ,111) 

CALl  VVALC W2V,W2,UW2V,BW2,111 ) 

C         NOPMALI7E  CONVOY  TRACK  FOR  EACH  UNIFORMLY 

C         DISTRIBUTED  SUBMARINE  STARTING  POSITION 

C  SAMPLE  SIZE  FOR  NORMAL  IS  DETERMINE  BY  THIS  LOOP 

DP  7  JJ=1.100D 
YP=NOPM(2C. )-(PL/2. ) 

C         DETERMINE  IF  DETECTION  OCCURS 

C         EACH  CPUNT  OF  THE  LOOP  INDEX  REPRESENTS  ONE  SAMPl E 

DP  P  J=1,JKL 

C  DETEPMIMF  S1J«  POSIT  (UNIFORMLY  DISTRIBUTED) 

C  INPUTS 

C  XP:  X-DISTANCE  TO  DATROL  LINE 

C  YP:  Y-COORDINATE  OF  LOWER  LIMIT  CF  PATRPL 

C  LINE 

C  PL:  LENGTH  OF  PATROL  LINE 

r  OUTPUTS 

f  PC:   RANGE  TO  SUB'S  INITIAL  POSIT 

C  BC:   BEARING  OF  SUB'S  INITIAL  POSIT 

r  YP  (VIA  COMMON/ONE/):   INITIAL  Y-COORDINATE 

r  sub 

r  P!V:  INITIAL  SUB  POSITION  VECTOR 

CALL  WYFP(XP,YPfPL,RO,BD,RlV) 
f  TTMF  TO  END  OF  FIRST  (RANDOM)  LFG 

IF(SGN)7CC,10«5,702 

C         TIM^  TC  GO  FCR  SUB  TO  REACH  LIMIT  OF  PATROL  LINE 
C  IF  TNITTAL  SUB  VELOCITY  IS  "DDWN" 
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7f\C\ 


702 
-"3 


Lr'3 


1.  0 1 


10? 


!  06 


11 


T]=( YP-YF ) /V 

or  t p  ? c ? 

TJMT  TT  00  FOR  SUP  TO  REACH  UPPFR  LIMIT  OF  PATPOl 

i  n,rr 

Tl=( YP+PL-YR)/V 
T  =  T1 

INITIAL T7F  TOTAL  TI^F 

TN=o. 

INPTT^Tf  THAT  THE  FIRST  (RANDOM!  L^G  IS  TO  BE 
F  XAMIMFP 

COMPUTE  LENGTH  OF  RANDOM  LEG 
GP  Tp  !C^ 

I  PGP  TP  INCPFMENT  LEO  CGUNT 
pr  ic     TM=1,MM 

CHECK  Tp  OETFR^INE  INITIAL  DIRFCTION  OF  SUB 

IF(SON)lfl ,10?, 102 

INITIAL  HTRFCTION  IS  UP:  USE  PARAMETERS  FPR  THE 
UP  LFG 

TW=TW? 

f/LL  CHPI (PWV,UW2V) 

w=  w? 

CALl  Cnpi(UV,W2V) 
C  £  LL  CPPI(VV,V2V) 
OC  TP  \rf 

INITIAL  DIRECTION  IS  GOWN:  USE  PARAMETERS  PnR  THF 
DTWN  LEG 

call  cnPnuvivf uwiv) 

TW=tvji 

W=W1 

C£LL     CnPI(  V>'V,W]  V) 

CALL     C°CI(VV,V1V) 

REVEPSE    DIRECTION    GF    NEXT    LEG 
SGN=-SON 

CCMOUTE    RELATIVE    LFNGTH    Oc    THIS    LEO 
FAIL     VMIIt  T(')WV,TFMV,TW) 

PtTFRMIME  VFFTGR  FOR  FINAL  POSITION  CM  THIS  LPG 

PALI     VADP(PlV,Tf:MVtPNVt  1  ) 

C£U     CPACHK(niV,PNVtlJWV,PMIN,INt  J) 

IMCPFME^T    LFG    COUNT 

SI EO=SLFG+l. 

DETFPMINE  I p  DETECTION  OCCURRED  ON  THIS  LEG 

IF(PMTM-R )3,3,11 

TN=TN+T 

IF(tn-tM^x)6,8,P 
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ft    IF(MM-1)1C5,9, 10 
10    CfNTINUF 
or   TO    8 
9    T=PL/V 
mm  =  i  OOCT 
IOA    TW1=T*W1 
tw?=T*W2 

no  Tr   ic^ 
?  nrT=rET-n. 

«    CPMTTMUF 

7  crNTlNiiF 

XJKL=JKl*K  00 
P=PFT/XJKL 
Xirn=SLFG/XJKL 

WFITF(6,6C1 )XLFG,U,RHC,Wl,W2,ALF,P 
'01    U=U+10. 

VFITF(f ,f C2I 
AO?  p^PM/lK/) 
TCn  ALF=ALF+PI/36, 

«,0?  FPPM£T( 1H1/////T16, 'PROBABILITY  OF  RANDOM  ENCOUNTER',/ 
l//T6,'DF  -I  MI TF  RANGE  L AW • /T8 , 'DETFCTI ON  RANGF:  »,F5.1, 
?*  MUFf  '//Tft,  'LENGTH  CF  SUB  &ATPOL  LINF:  •fF5.1,«  mil' 
?,  »ES*//T6 , 'X-DISTANCF  TO  PATROL  LI NE: • , F5 . 1 , •  MILFS'// 
tit  ,  "SUBMARINE  SPEED: • fF5. 1. •  KNOTS  '  //T6,  •  NO.  0<=  SAMPL' 
C,'PS  TAKEN  C0P  EACH  SCENARIO:  • ,I4///T3f«  AVG.  NO.  SH' 
*,'IP   FPFFT  ■ 

^♦•PfLATIVF  SPFFD   ALFA    PRHB.  DETECT  ION • /T5 ,' L FGS    • 
p,»  SPECP  P^tio  UPLFG   DCWNL5GV/1 
•SOI  FrRM£T(T^fF5.2t?X,3(lXfF5.l  > ,  2X  ,  F5.  1 ,  3X,  F6.  4,  7X,  F6.  4  ) 
105  STOP 
END 


SUBROUTINE  WYEB ( XP, YP , PL ,R0 , BO, A) 

C  INPUTS 

C  YB=Y-POSITION  OF  SUB  ON  PATROL  LINF 

C  YB  IS  UNIFORMLY  DISTRIBUTED 

C  SGN*INIT!AL  DIRECTION  (EQUALLY  LIKELY)  OF  SUB 

C  VFLOCITY 

CCWMON/ONE/SGNfYB 
DIMENSION  A<2)  ,B<2) 

C         DETFPM^C  SUB  POSIT  RANDOMLY 

YP=YP+UPN( 1 )*PL 

C         DETERMINE  INITIAL  DIRECTION 

*  IF(UPN(l)-r.5l5,6,7 

a  or  to  ^ 

C         UP 

*  scm=i. 

GO  TO  o 

C  DOWN 

7  SGN*-1, 

r         DETERMINE  INITIAL  RANGF  AND  BEARING  TH  SUB 

3  M 1)=XP 
A(2)=YB 
CALL  WAL(  A,R0,P,B0,101) 
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RFTURN 
END 


SUBROUTINE  CPACHK ( Rl, RN,UWV,RCP A, IN, J ) 

DIMENSION  Rim,RN<l),UWV<l),X(2),Y(21,CPA<2l,Z<2) 
COMMON/TWO/Y,N 

C         R1=REFERENCE  VECTOR 

C         RN=TERMINAL  VECTOR 

C         UWV=UNIT  VECTOR  OF  DRM 

C         COMPUTE  DOT  PRODUCT:  R1*UWV=X 
CALL  VADD(R1,UWV,X,0) 


IF(XU))1,2,2 


C         CPA  WILL  OCCUR  ON  THIS  LEG:  COMPUTE  THE  DOT 
C         PRODUCT:  RN*UWV=X 

1  CALL  VADD<RN,UWV,X,0) 
IF(X(1))3,3,4 

C         CPA  WILL  NOT  OCCUR  ON  THIS  LEG:  CPA*R1   COMPUTE 
C         CPA 

2  CALL  VVAL(P1,RCPA,X,D,100) 
RFTURN 

C         CPA  WAS  NOT  REACHED  ON  THIS  LEG:   CPA*RN:  COMPUTE 
C         CPA 

3  CALL  VVAL(RN,RCPA,X,D,100) 
RFTURN 

C         DETERMINE  CPA:  ASSUME  LEG  HAS  LENGTH<  250  MILES 

4  DO  5  1=1,500 

C         USE  HALF-MILE  INCREMENTS 

P=I/2« 
C         INCREMENT  DRM  BY  1/2  MILE 

CALL  VMULT(UWV,Y,P) 
C         ADD  TO  REFERENCE  VECTOR 

CALL  VADD(R1,Y,X,1) 
C         COMPUTE  DOT  PRODUCT   X  DOT  UWV=Z 

CALL  VADD(X,UWV,Z,0) 

C         IF  Z  CHANGES  SIGN.  OR  BECOMES  ZERO  VECTOR  X=CPA: 
C         OTHERWISE  TRY  AGAIN 

IF(Z(1M  5,6,6 

5  CONTINUE 

C         DRM  LEG>250  MILES 

WPITE(6,60)IN,  J 
60  FOPMAT(//10X,«LEG»,1X,I4,1X,»OF  SAMPLE* , 15 , IX, • IS 
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1GFFATFP  THAN  250  HILFSV/I 
RF^UPN 
6  CALL  VVAL(X,RCPA,Y,D,100) 

PFTURN 


SLPPFUTINF  VMULM  4fB,X) 

C  THIS  SUBROUTINE  MULTIPLIES  THE  VECTGP  A  BY  THF 

C         SCALAR  X  TO  FORM  THF  VECTOR  B 

niMFNSTCN  A< I) ,B( II 
00  1  1=1,2 
1  B( I)=X*A(I I 

BFTUPN 

TMO 


SUBPOUTINF  VVAL<A,VL,UV,ANGLF,M) 

C  THIS  SUBROUTINE  DETERMINES  THE  FOLLOWING 

C  VALUES  OF  A: 

C  A=INPUT  VECTOR  (SUPPLIED  BY  THE  CALLING  PROGRAM) 

C  VL*MAGNITUDE  OF  A 

C  UV*UNIT  VECTOR  PARALLEL  TO  A 

C  ANGLF*ANGLE  OF  A  MEASURED  CCW  FROM  X-AXIS 

DIMENSION  A(l> ,UV<1) 
IF(M-100)2,3,4 
<*    IMM-11015,6,7 

C         COMPUTE  VECTOR  LENGTH  ONLY 

3  VL*SQPT(A<1I**2^A(2)**2I 

RETURN 

C         COMPUTE  ANGLE  OF  VECTOR  ONLY 

2  ANGLE*ATAN2<A(2),A(l)) 

RFTURN 

C         COMPUTE  BOTH  LENGTH  AND  ANGLE  OF  VECTOR 

•5  Vl=SQRT(A(l)**2*A(2)**2) 

GO  TO  2 

C         GENERATF  UNIT  VECTOR  ONLY 

6  VL*SQRT(A(1)**2*A(2)**2) 
DO  1  1*1.2 

1  UV(I )=A(f l/VL 
RETURN 

C         COMPUTE  ANGLE  AND  GENERATE  UNIT  VECTOR 

7  ANGLE=ATAN2(A(2)  ,A(D) 
GO  TO  6 

END 
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SUBROUTINE  VINS(A,X,Y) 

C         THIS  SUBROUTINE  INSERTS  THE  X  AND  V  VECTOR 
C         COMPONENTS  INTO  ARRAY  A 

DIMENSION  A(l) 

A(1)=X 

A(2)=Y 

RETURN 

END 


SUBROUTINE  COPI(AfB) 

C         THIS  SUBROUTINE  GENERATES  THE  VECTOR  A  SUCH  THAT 
C         FOR  A  GIVEN  VECTOR  B,  A-B 

DIMENSION  A(1)VB(1I 
DO  1  1-1.2 
I  A(I)-B(T1 

RETURN 
END 


SUBROUTINE  VAOD( A ,B,R ,N) 

C  THIS  SUBROUTINE  DOES  THE  FOLLOWING  VECTOR 

C  OPERATIONS: 

C  IF  N>0:  R-A+8 

C  IF  N<0:  R=A-B 

C  IF  N=0:  THE  DOT  PRODUCT  OF  A  AND  B  IS  FORMED  AND 

C  RESULT  IS  PLACED  IN  THE  FIRST  MEMBER  OF  ARRAY  R 

DIMENSION  A(l)  ,B(1)  ,R(1) 
IF(N)1,2,3 

C         ADD  VECTORS 

3  DO  4  1=1,2 

A  M1)-A<U*B<II 

RETURN 

C         SUBTRACT  VECTORS 

1  DO  5  1-1.2 

5  R(I)-A(I)-B(I1 

RETURN 

C         COMPUTE  DOT  PRODUCT  OF  VECTORS 

2  DO  6  1-1.2 

6  R(II-A(I)*B(I) 
R(1)=R(1)+P(2) 

PEniRN 
END 
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THlr>  SlJPRr)UTTN"  GFNE"ATES  MOfSIGWA  SQUARED) 
PANDPM  VARIABLES  *Y  THF  ROX-MUtLLFP  METHTO 

TWC  INPF°FNPENT  NORMALLY  DISTRIBUTED  RANDO* 
VAPURlcS  APE  PRODUCED  EACH  TlMr  THIS  SUBROUTINE 
T  S  CMin 


P  =  STANCARO  CF VI  AT  I  ON 


PFU  FUNCTION 


D^»  K.TWPFI/lf* 


NORM (0 J 

23318  53/ 
CC  TP  ( 1-21  ,N 
p  =  SOPT(-?.*ALOC(URN(n  »  ) 
APG=TWPPI*(URM( 1) I 

Nrpp=p*cr.c  (  arc  )*n 

X?=p*cim( APG)*P 

N  =  ? 

Rf  TURN 

MfiPM-X? 

M*l 

PFTURN 

EN'P 
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APPENDIX  E 
A.   GEOMETRIC  INTEGRATION 

Figure  7  is  a  graph  of  a  family  of  curves  plotted  from 
Table  E-la  which  represents  the  conditional  probability 
that  detection  occurs  as  a  function  of  the  crossing  point 
displacement  from  the  patrol  line  center.   The  normal 
density  function,  with  mean  zero  and  standard  deviation  20 
miles,  is  also  plotted  in  the  figure.   The  upper  set  of 
numbers  along  the  abscissa  represents  the  normalized  values, 
i.e.  , 

x'  =  x/20 

where  x'  is  the  normalized  value  of  x.   The  product  of  the 
normal  density  curve  with  each  of  the  conditional  proba- 
bility curves  is  tabulated  in  Table  E-lb,  and  plotted  as 
a  family  of  curves  in  Figure  9.   If  a  grid  of  0.1  milli- 
meter squares  is  superimposed  on  this  figure,  the  area 
represents  one-half  of  the  integral  of  the  function.   The 
area  under  each  curve  was  determined  by  the  grid-counting 
method,  and  was  also  determined  by  use  of  the  Amsler  Polar 
Planimeter.   Each  curve  was  traced  five  times  with  the 
planimeter  stylus,  and  an  average  planimeter  dial  reading 
taken.   The  results  of  these  two  methods  of  geometric 
integration  are  tabulated  in  Table  9.   The  average  plani- 
meter dial  readings  were  normalized  to  0.998  (the  value  of 
P(D)  for  p  =  0)  using  the  following  equation: 
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TABLE  E-l 

CONDITIONAL  PROBABILITY  OF  DETECTION 
GIVEN  CROSSING  POINT  IS  DISPLACED 
Y  MILES  FROM  PATROL  CENTER 


a.   NO  ASSUMPTION  ON  DISTRIBUTION  OF  CROSSING  POINTS 


P(D|Y 

IARINE^S 

=  y) 

y 

SHIP-SUBI* 

IPEED  RATIO  (p) 

p  =  0 

P-l 

p  =  1.5 

P  =  2 

P  =  3 

p=5 

p=9 

P=oo 

0 

1 

1 

.9354 

.8000 

.6250 

.5104 

.5021 

.5000 

5 

1 

1 

.9260 

.7844 

.6250 

.5260 

.5021 

.5000 

10 

1 

1 

.8687 

.7635 

.6177 

.5469 

.5021 

.5000 

15 

1 

.9781 

.8062 

.7125 

.6104 

.5604 

.5292 

.5000 

20 

1 

.9365 

.7437 

.6500 

.5969 

.5604 

.5292 

.5000 

25 

1 

.8844 

.6812 

.5875 

.5552 

.5469 

.5292 

.5000 

30 

1 

.8219 

.6187 

.5302 

.5135 

.5052 

.5010 

.5000 

50 

1 

.5719 

.3844 

.3635 

.3469 

.3385 

.3344 

.3333 

70 

1 

.3219 

.2177 

.1969 

.1802 

.1719 

.1677 

.1667 

b.   NORMAL  DISTRIBUTION  (STANDARD  DEVIATION  20  MILES) 
OF  CROSSING  POINTS  ASSUMED 


p(d|y  = 

y) 

y 

<My)* 

SHII 

P 

»  SUBMARINE  SPEED  RATIO  ( p) 

p  =  0 

P-l 

p=1.5 

P*2 

P=3 

P=5 

p=9 

P  =  oo 

0 

.3989 

.3989 

.3731 

.3191 

.2493 

.2036 

.2003 

.1994 

5 

.3867 

.3867 

.3581 

.3033 

.2417 

.2034 

.1942 

.1934 

10 

.3521 

.3521 

.3059 

.2688 

.2175 

.1926 

.1768 

.1760 

15 

.3011 

.2945 

.2428 

.2145 

.1838 

.1687 

.1574 

.1506 

20 

.2420 

.2266 

.1800 

.1573 

.1444 

.1356 

.1281 

.1210 

25 

.1826 

.1615 

.1244 

.1073 

.1014 

.0999 

.0966 

.0913 

30 

.1295 

.1064 

.0812 

.0687 

.0665 

.0654 

.0649 

.0648 

50 

.0175 

.0100 

.0067 

.0064 

.0061 

.0059 

.0058 

.0058 

70 

.0009 

.0006 

.0002 

.0002 

.0004 

.0002 

.0002 

.0003 

♦Extracted  from  CRC  Standard  Mathematical  Tables, 
1968,  pp.  549-555. 


16th  ed. 
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Area  under  Curve  =  (0.9980  x  average  reading) /762    (2) 

where  the  constant  762  is  the  average  planimeter  dial 
reading  obtained  by  tracing  the  curve  for^ p  =  0. 
The  values  contained  in  Table  E-2  for  the  column  under  the 
heading  "Planimeter  P(D)"  were  computed  by  equation  (2) 
using  a  Friden  electronic  calculator. 

TABLE  E-2 


RESULTS  OF  TWO  GEOMETRICAL 
INTEGRATION  METHODS 


SPEED  RATIO 
<<  P). 


P  (DETECTION) 


GRID  COUNTING       AMSLER  POLAR 

PLANIMETER  * 


0  0.9980  0.9980 

1  0.9415  0.9319 
1.5  0.7895  0.7922 

2  0.6816  0.6828 

3  0.5781  0.5797 
5  0.5233  0.5234 

9  0.5008  0.5011 

0.4885  0.4883 


The  planimeter  dial  readings  were  normalized  so  that 
762  increments  on  the  dial  represented  the  known  area  0.99  80, 
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APPENDIX  F 


LIMITING  CASE  FOR  P (DETECTION) 


C 

C 

C  LIMIT 

C 

C 

J=0 

R=30. 

PL=120. 

YP=40. 

WRITE (6, 600 )R, PL 

Z=YP+PL 

DO  1  1=1,41 

IF(YP-R) 3,2,2 

2  P=0 

GO  TO  4 

3  IF (YP+R) 5,5,6 

5  P=(2.*R)/PL 
GO  TO  4 

6  Y=AMIN1(ABS (YP) ,ABS(Z) ) 
S=(YP*Z) 

W=S/ABS (S) 
P=(R-W*Y)/PL 

4  J=J+1 
IF(J-6) 7,8,9 

8  J=0 
WRITE(6,602) 

7  WRITE(6,601)YP,P 
1  YP=YP-2.5 

9  STOP 

600  FORMAT (1H1,1 5 (/)  ,20X, 'LIMITING  CASE  FOR  P (DETECTION) '/ 
1//5X, 'DETECTION  RADIUS  (DEFINITE  RANGE  LAW):  ',F4.0//5 
2X,  'PATROL  LENGTH:  ' ,F4 . 0///7X, ' YP ' , 8X, 'P (D) ' //) 

601  FORMAT (5X,F5 . 1 , 5X,F6 . 4) 

602  FORMAT  (/) 
END 
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TABLE  F-l 
LIMITING  CASE  FOR  P  (DETECTION) 

DETECTION  RADIUS  (DEFINITE  RANGE  LAW):  30. 
PATROL  LENGTH:   120. 

YP  P(D) 

40.0  0.0 

37.5  0.0 

35.0  0.0 

32.5  0.0 

30.0  0.0 

27.5  0.0208 

25.0  0.0417 

22.5  0.0625 

20.0  0.0833 

17.5  0.1042 

15.0  0.1250 

12.5  0.1458 

10.0  0.1667 

7.5  0.1875 

5.0  0.2083 

2.5  0.2292 

0.0  0.2500 

-2.5  0.2708 

-5.0  0.2917 

-7.5  0.3125 

-10.0  0.3333 

-12.5  0.3542 

-15.0  0.3750 

-17.5  0.3958 

-20.0  0.4167 

-22.5  0.4375 

-25.0  0.4583 

-27.5  0.4792 

-30.0  0.5000 

-32.5  0.5000 

-35.0  0.5000 

-37.5  0.5000 

-40.0  0.5000 

-42.5  0.5000 

-45.0  0.5000 

-47.5  0.5000 

-50.0  0.5000 

-52.5  0.5000 

-55.0  0.5000 

-57.5  0.5000 

-60.0  0.5000 
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Ship-Submarine  Speed  Ratio:  °° 

Submarine  Patrol  Line  Length: 120 

miles 
Detection  Range  (Definite  Range 
Law) :   30  miles 


0     10    20     30    40     50    60     70    80     90 
Crossing  Position  in  Miles  (from  Patrol  Line  Center) 


Figure  F-l.   Probability  of  Detection  vs.  Crossing 
Position  (Perpendicular  Tracks)  for  the 
Limiting  Case 
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Attacks  (U) ,  p.  31,  February  28,  1955.  (C) 

2.  Operations  Evaluation  Group,  Office  of  the  Chief  of 
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4.  Operations  Evaluation  Group,  Office  of  the  Chief  of 

Naval  Operations,  Report  No.  56,  Search  and  Screening, 
by  B.  0.  Koopman,  p.  106-107,  1946. 

5.  Ibid,  p.  107,  Figure  8. 
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p.  279,  McGraw-Hill,  1956. 
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